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Abstract 

Relativistic energy density functionals (EDF) have become a standard tool 
for nuclear structure calculations, providing a complete and accurate, global 
description of nuclear ground states and collective excitations. Guided by 
the medium dependence of the microscopic nucleon self-energies in nuclear 
matter, semi-empirical functionals have been adjusted to the nuclear mat- 
ter equation of state and to bulk properties of finite nuclei, and applied to 
studies of arbitrarily heavy nuclei, exotic nuclei far from stability, and even 
systems at the nucleon drip-lines. REDF-based structure models have also 
been developed that go beyond the static mean-field approximation, and 
include collective correlations related to the restoration of broken symme- 
tries and to fluctuations of collective variables. These models are employed 
in analyses of structure phenomena related to shell evolution, including de- 
tailed predictions of excitation spectra and electromagnetic transition rates. 

Keywords: Nuclear structure, Nuclear Density Functional Theory and 
extensions. Collective correlations 



1. Introduction 

Nuclear energy density functionals (EDF) presently provide the most 
complete and accurate description of ground-state properties and collective 
excitations over the whole nuclide chart. Among the microscopic approaches 
to the nuclear many-body problem, probably no other method achieves com- 
parable global accuracy at the same computational cost, and it is the only 

Preprint submitted to Progress m Particle and Nuclear Physics January 12, 2013 



one that can describe the evolution of structure phenomena from relatively 
light systems to superheavy nuclei, and from the valley of /3-stability to the 
particle drip-lines [l|, [3] . 

In practical implementations the EDF framework is realized on two spe- 
cific levels. The basic implementation is in terms of self-consistent mean-field 
(SCMF) models, in which an EDF is constructed as a functional of one-body 
nucleon density matrices that correspond to a single product state ~ Slater 
determinant of single-particle or single-quasiparticle states. The SCMF ap- 
proach to nuclear structure is analogous to Kohn-Sham density functional 
theory (DFT) 0, fl. DFT enables a description of quantum many-body 
systems in terms of a universal energy density functional. Universal in the 
sense that, for a given inter-particle interaction, it has the same functional 
form for all systems. Nuclear SCMF models effectively map the many-body 
problem onto a one-body problem, and the exact EDF is approximated by 
simple, mostly analytical, functionals of powers and gradients of ground- 
state nucleon densities and currents, representing distributions of matter, 
spins, momentum and kinetic energy. In principle the nuclear EDF can in- 
corporate short-range correlations related to the repulsive core of the inter- 
nucleon interaction, and long-range correlations mediated by nuclear reso- 
nance modes. Even though it originates in the effective interaction between 
nucleons, a generic density functional is not necessarily related to any given 
nucleon- nucleon (NN) potential and, in fact, some of the most successful 
modern functionals are entirely empirical. Of course it would be desirable 
to have a fully microscopic foundation for a universal density functional, 
and this is certainly one of the major challenges for the framework of nu- 
clear EDFs. Because it includes correlations, the self-consistent Kohn-Sham 
approach goes beyond the Hartree-Fock approximation, and it also has the 
advantage of being a local scheme. Its usefulness, however, crucially depends 
on our ability to construct accurate approximations for the most important 
art of the functional, that is, the universal exchange-correlation functional 
In an "ab initio" approach one might start from a Hamiltonian that 
describes two-nucleon and few-body scattering and bound-state observables 
0] , or an effective field theory of low-energy in-medium NN interactions can 
be used to build approximations to the exact exchange-correlation functional 
However, even if a fully microscopic EDF is eventually developed, the 
parameters of that functional will still have to be fine tuned to structure 
data of finite nuclei. This is because data on nucleon-nucleon scattering 
and few-nucleon systems, or gross properties of infinite nuclear matter, can- 
not determine the density functional to a level of accuracy necessary for a 
quantitative description of medium-heavy and heavy nuclei. 
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When considering applications, however, an important challenge for the 
framework of EDF is the systematic treatment of collective correlations re- 
lated to restoration of broken symmetries and fluctuations in collective co- 
ordinates. A static nuclear EDF is of course characterized by symmetry 
breaking ~ translational, rotational, particle number, and can only provide 
an approximate description of bulk ground-state properties. To calculate 
excitation spectra and electromagnetic transition rates in individual nuclei, 
it is necessary to extend the Kohn-Sham EDF framework, that is the SCMF 
scheme, to include correlations that arise from symmetry restoration and 
fluctuations around the mean-field minimum. Collective correlations are 
sensitive to shell effects, display pronounced variations with particle number 
and, therefore, cannot be incorporated in a universal EDF. On the second 
level that takes into account collective correlations through the restoration of 
broken symmetries and configuration mixing of symmetry-breaking product 
states, the many-body energy takes the form of a functional of all transition 
density matrices that can be constructed from the chosen set of product 
states. This set is chosen to restore symmetries or/and to perform a mix- 
ing of configurations that correspond to specific collective modes using, for 
instance, the (quasiparticle) random-phase approximation (QRPA) or the 
Generator Coordinate Method (GCM). The latter includes correlations re- 
lated to finite-size fluctuations in a collective degree of freedom, and can be 
also used to restore selection rules that are crucial for spectroscopic observ- 
ables. 

An important class of nuclear structure models belongs to the frame- 
work of relativistic energy density functionals. In particular, a number of 
very successful relativistic mean-field (RMF) models have been constructed 
based on the framework of quantum hadrodynamics (QHD) [^. There 
are important advantages in using functionals with manifest covariance [lol |. 
The most obvious is the natural inclusion of the nucleon spin degree of free- 
dom, and the resulting nuclear spin-orbit potential which emerges automat- 
ically with the empirical strength in a covariant formulation. The consistent 
treatment of large, isoscalar, Lorentz scalar and vector self-energies provides 
a unique parametrization of time-odd components of the nuclear mean-field, 
i.e. nucleon currents, which is absent in the non-relativistic representation 
of the energy density functional. The empirical pseudospin symmetry in nu- 
clear spectroscopy finds a natural explanation in terms of relativistic mean 
fields On a microscopic level, it has been argued that a covariant 
formulation of nuclear dynamics manifests the true energy scales of QCD 
in nuclei, and is consistent with the nonlinear realization of chiral symme- 
try through the implicit inclusion of pion-nucleon dynamics in the effective 
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nucleoli self-energies. A covariant treatment of nuclear matter provides a 
distinction between scalar and four- vector nucleon self energies, leading to 
a very natural saturation mechanism. 

RMF-based models have been very successfully employed in analyses 
of a variety of nuclear structure phenomena, not only in nuclei along the 
valley of /3-stability, but also in exotic nuclei with extreme isospin values 
and close to the particle drip lines. Applications have reached a level of so- 
phistication and accuracy comparable to the non-relativistic Hartree-Fock- 
Bogoliubov approach based on Skyrme functionals or Gogny effective inter- 
actions 0, 12. lid. 14 1 . Relativistic EDFs have mostly been applied in the 
description of ground-state properties and excitation energies of giant reso- 
nance at the self-consistent mean-field level, taking into account pairing cor- 
relations in open-shell nuclei in the Hartree-(Fock)-Bogoliubov framework, 
and performing consistent (Q)RPA calculations of small-amplitude collective 
motion. However, for relativistic structure models to make detailed spectro- 
scopic predictions, symmetries broken by the static nuclear mean field must 
be restored, and fluctuations around the mean-field minimum must be taken 
into account. While symmetry restoration and configuration mixing calcu- 
lations have routinely been applied with non-relativistic density functionals 
since many years, it is only more recently that this type of structure models 
have been developed using relativistic density functionals jl5l.ll6l.ll7l.ll8l.ll9l|. 

In this work we review recent advances in the framework of relativistic 
EDFs and, in particular, the latest extensions that include the treatment 
of collective correlations. Section [2] introduces the general framework of 
REDFs. In Section [3] we review a class of semi-empirical functionals and, 
in particular, the functional DD-PCl that will be used in illustrative cal- 
culations throughout this work. For a quantitative analysis of open-shell 
nuclei it is necessary to consider pairing correlations, and in Section U] we 
introduce the relativistic Hartree-Bogoliubov model for triaxial nuclei, with 
a separable pairing interaction. The treatment of collective correlations is 
reviewed and illustrated with a number of examples in Section [5] (symmetry 
restoration and configuration mixing calculations) and in Section [6] (collec- 
tive Hamiltonian in five dimensions). Section [7] summarizes the results and 
ends with an outlook for future studies. 



2. Relativistic Energy Density Functionals 

2.1. Quantum Hadrodynamics 

In conventional quantum hadrodynamics (QHD) jl, 0] a nucleus is de- 
scribed as a system of Dirac nucleons coupled to exchange mesons through 
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an effective Lagrangian. Tlie isoscalar scalar a meson, tlie isoscalar vector 
u meson, and tlie isovector vector p meson build the minimal set of me- 
son fields that, together with the electromagnetic field, is necessary for a 
description of bulk and single-particle nuclear properties. In the mean-field 
approximation the meson-field operators are replaced by their expectation 
values in the nuclear ground state. In addition, a quantitative treatment of 
nuclear matter and finite nuclei necessitates a medium dependence of effec- 
tive mean-field interactions that takes into account higher-order many-body 
effects. A medium dependence can either be introduced by including non- 
linear meson self-interaction terms in the Lagrangian, or by assuming an 
explicit density dependence for the meson- nucleon couplings. The former 
approach has been adopted in the construction of several successful phe- 
nomenological RMF interactions, for instance, the very popular NL3 [201 ] . 
or the more recent PKl, PKIR 21] and FSUGold [i^] parametrizations of 
the effective Lagrangian. In the latter case, the density dependence of the 
meson-nucleon vertex functions can be parameterized starting from micro- 
scopic Dirac-Brueckner calculations of symmetric and as ym metric nuclear 



matter [23, 1 24 . |25| or it can be fully phenomenological 26|, 1271 . 1281 ]. with 



parameters adjusted to data on finite nuclei and empirical properties of 
symmetric and asymmetric nuclear matter. 

At the energy scale characteristic for nuclear binding and low-lying ex- 
cited states, meson exchange (cr, u), p, . . .) is just a convenient represen- 
tation of the effective nuclear interaction. The exchange of heavy mesons 
is associated with short-distance dynamics that cannot be resolved at low 
energies, and therefore in each channel (scalar- isoscalar, vector- isoscalar, 
scalar-isovector, and vector-isovector) meson exchange can be replaced by 
the corresponding local four-point (contact) interactions between nucleons. 
The self-consistent relativistic mean-field framework can be formulated in 
terms of point-coupling nucleon interactions. When applied in the descrip- 
tion of finite nuclei, relativistic mean-field point-coupling (RMF-PC) models 
l33| produce results that are equivalent to those obtained in the 



meson exchange picture. Of course, also in the case of contact interactions, 
medium effects can be taken into account by the inclusion of higher-order 
interaction terms, for instance, six- nucleon vertices {'tptp)^, and eight-nucleon 
vertices {tptl^)^ and [('i/'7^t/')(V'7'^V')]^5 or it can be encoded in the effective 
couplings, i.e. in the density dependence of strength parameters of the in- 
teraction in the isoscalar and isovector channels. Although a number of 
point-coupling models have been developed over the years, it is only more 
recently that phenomenological parametrizations have been adjusted and 
applied in the description of finite nuclei on a level of accuracy comparable 
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to that of standard meson-exchange effective interactions [33|, \34 

The relation between the two representations: finite-range (meson ex- 
change) and zero-range (point-couphng), is straightforward in nuclear mat- 
ter because of constant nucleon scalar and vector densities. The Klein- 
Gordon equations of the meson-exchange model with meson masses and 
density-dependent couplings ff<^(/o), are replaced by the corresponding point- 
coupling interaction terms with strength parameters g^/fn^- In finite nuclei, 
however, the problem is not so simple. Because of the radial dependence 
of the densities, the expansion of the meson propagator in terms of ^/fn^ 
leads to an infinite series of gradient terms. In practice this series has to 
be replaced by a finite number of terms with additional phenomenological 
parameters adjusted to low-energy data. A number of studies have shown 
that, both for finite-range and for point-coupling mean-field models, the em- 
pirical data set of ground-state properties of finite nuclei can determine only 
a relatively small number of parameters in the general expansion of the effec- 
tive Lagrangian in powers of the fields and their derivatives. It is therefore 
not a priori clear how to select the set of point-coupling interaction terms 
that will describe structure properties at the same level of accuracy as the 
meson-exchange models. The mapping of a phenomenological finite-range 
interaction with density-dependent meson-nucleon couplings (DD-ME2) on 
the zero-range (point-coupling) relativistic mean-field framework was con- 
sidered in Ref. [35|]. A family of point-coupling effective interactions was 
constructed with different values of the strength parameter of the isoscalar- 
scalar derivative term. In the meson-exchange picture this corresponds to 
different values of the cr-meson mass. The parameters of the isoscalar-scalar 
and isovector-vector channels of the point-coupling interactions were ad- 
justed to nuclear matter and ground-state properties of finite nuclei. By 
comparing results for infinite and semi- infinite nuclear matter, ground-state 
masses, charge radii, and collective excitations, constraints were placed on 
the parameters of phenomenological point-coupling relativistic effective in- 
teraction. 



2.2. Elements of Density Functional Theory 

Density Functional Theory (DFT) is one of the most popular and success- 
ful "ab initio" approaches to the structure of quantum many-body systems 
(atoms, molecules, solids). Probably no other method achieves compara- 
ble accuracy at the same computational cost. The basic concept is that 
the ground-state properties of a stationary many-body system can be repre- 
sented in terms of the ground-state density alone. Since the density p(r) is 
a function of only three spatial coordinates, rather than the 3N coordinates 
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of the N-body wave function, DFT is computationally feasible even for large 
systems. 

Most practical applications of Density Functional Theory use the effec- 
tive single-particle Kohn-Sham (KS) equations 0,0], introduced for an aux- 
iliary system of N non-interacting particles. According to the Hohenberg- 
Kohn theorem [s^, there exists a unique energy functional 



Es[p\=TM + j d\vs{v)p{v), (1) 



for which the variational equation yields the exact ground-state density 
/9s(r). Ts[p\ is the universal kinetic energy functional of the non-interacting 
system. The KS scheme is based on the following assertion: for any in- 
teracting system, there exists a unique local single-particle potential t's(r), 
such that the exact ground-state density of the interacting system equals 
the ground-state density of the auxiliary non-interacting system: 

N 



p{t) = p,(r) =^\(t)i{T) 



expressed in terms of the N lowest occupied single-particle orbitals - solu- 
tions of the Kohn-Sham equations: 

[-VV2m + Vs{y)\ 0,(r) = ^^^.(r) . (3) 

The uniqueness of Vs{v) follows from the Hohenberg-Kohn theorem and 
the single-particle orbitals are unique functionals of the density: </>i(r) = 
</'i(M;r). 

For a self-bound system like the atomic nucleus, the energy functional 
can be decomposed into three separate terms: 

F[p] = TM + Eh[p\ + EAp\. (4) 

where Tg is the kinetic energy of the non-interacting A-nucleon system. Eh 
is a Hartree energy, and E^c denotes the exchange-correlation energy which, 
by definition, contains everything else - all the many-body effects. The 
corresponding local exchange-correlation potential is defined by: 

and thus 

Vs [p\ (r) = VH [p] (r) + Vx-c [p] (r ) . (6) 
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Since the effective potential depends on tlie ground-state density, tfie system 
of equations (H]), ([3l), and ^ has to be solved self-consistently. This is the 
Kohn-Sham scheme of density functional theory [3] . By including correlation 
effects the KS framework goes beyond the Hartree-Fock approximation but, 
in addition, it has the advantage of being a local scheme. It is clear, however, 
that the usefulness of the Kohn-Sham scheme crucially depends on our abil- 
ity to construct accurate approximations to the exact exchange-correlation 
energy. The true exchange-correlation energy functional is universal, i.e. 
given the inter-particle interaction, it has the same functional form for all 
systems. One possible approach is to develop E^c from first principles by 
incorporating known exact constraints. Another is empirical, a parametric 
ansatz is optimized by adjusting it to a set of data. 

The Hohenberg-Kohn theorem and the self-consistent Kohn-Sham scheme 
are straightforwardly extended to the relativistic domain . The relativistic 
Kohn-Sham equation for the auxiliary non-interacting system is represented 
by the single-particle Dirac equation with a local four-potential that depends 
on the ground-state four-current. 

2. 3. Empirical nuclear density functionals 

At the energy and momentum scales characteristic of nuclei, the only 
degrees of freedom that have to be taken into account explicitly in the de- 
scription of many-body dynamics are pions and nucleons. The behavior of 
the nucleon-nucleon (NN) interaction at long and intermediate distances is 
determined by one- and two-pion exchange processes. As already empha- 
sized, short-distance dynamics cannot be resolved at low energies that char- 
acterize nuclear binding and, therefore, it is represented by local four-point 
(contact) NN interactions, with low-energy (medium-dependent) parameters 
adjusted to nuclear data. These concepts of effective field theory and density 
functional theory methods have recently been used to derive a microscopic 
relativistic energy density functional framework constrained by in-medium 



QCD sum rules and chiral symmetry [37l . l38l |. The density dependence of 
the effective nucleon-nucleon couplings is determined from the long- and 
intermediate-range interactions generated by one- and two-pion exchange 
processes. They are computed using in-medium chiral perturbation the- 
ory, explicitly including A (1232) degrees of freedom [39]. Regularization 
dependent contributions to the energy density of nuclear matter, calculated 
at three-loop level, are absorbed in contact interactions with parameters 
representing unresolved short-distance dynamics. 

In this work we review a class of relativistic energy density functionals 



(REDFs) similar to that introduced in Refs. [37|, l38( but, instead of using 
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low-energy QCD constraints for the medium dependence of the parameters, 
a phenomenological ansatz is adjusted exclusively to data on nuclear ground 
states. This empirical approach, although guided by microscopic nucleon 
self-energies in nuclear matter, gives us more freedom to investigate in de- 
tail the relationship between global properties of a nuclear matter equation 
of state (volume, surface, and asymmetry energies) and the corresponding 
predictions for properties of finite nuclei. 

The basic building blocks of a relativistic nuclear energy density func- 
tional are the densities and currents bilinear in the Dirac spinor field ijj of 
the nucleon: 

tpOrTip , Or G {l,Ti} , r G {l,-f^J.,J5,75l^L,o■^J.u} ■ (7) 

Here are the isospin Pauli matrices and F generically denotes the Dirac 
matrices. The nuclear ground-state density and energy are determined by 
the self-consistent solution of relativistic linear single-nucleon Kohn-Sham 
equations. To derive those equations it is useful to construct an interac- 
tion Lagrangian with four-fermion (contact) interaction terms in the various 
isospace-space channels: 

isoscalar-scalar: (V'V')'^ 

isoscalar- vector: (■(/'7^-(/;)(V'7^V') 

isovector-scalar: {'ipTip) ■ (tpTip) 

isovector-vector: [ipT^^ijj) ■ (ipfj^^l;) . 

Vectors in isospin space are denoted by arrows. A general Lagrangian can 
be written as a power series in the currents ipOrTip and their derivatives, 
with hi ghe r-order terms representing in-medium many-body correlations 



[29l . l3d . l31l . |32| . l33l |. The problem however, as already emphasized, is that 
the empirical data set of bulk and single-particle properties of finite nuclei 
can only constrain a relatively small set of parameters in the general ex- 
pansion of an effective Lagrangian. An alternative, that directly leads to 
linear single-nucleon Kohn-Sham equations, is to construct a Lagrangian 
with second-order interaction terms only, with many-body correlations en- 



coded in density-dependent coupling functions [37|, |38|. In complete anal- 
ogy to the successful meson-exchange RMF phenomenology, in which the 
isoscalar-scalar a meson, the isoscalar-vector cj meson, and the isovector- 
vector p meson build the minimal set of meson fields that is necessary for 
a quantitative description of nuclei, an effective Lagrangian that includes 
the isoscalar-scalar, isoscalar-vector and isovector-vector four-fermion inter- 
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actions reads: 

C = '(/5(z7 • d — m)ijj 

-^55(W)(W)-e^7-^^^^V'- (8) 

In addition to the free-nucleon Lagrangian and the point-coupHng inter- 
action terms, when apphed to nuclei, the model must include the cou- 
pling of the protons to the electromagnetic field. The derivative term in 
Eq. ([8|) accounts for leading effects of finite-range interactions that are cru- 
cial for a quantitative description of nuclear density distribution, e.g. nuclear 
radii. Similar interactions can be included in each space-isospace channel, 
but in practice data only constrain a single derivative term, for instance 
l's{dp^'^){d^^^)- The inclusion of an adjustable derivative term only in 
the isoscalar-scalar channel is consistent with conventional meson-exchange 
RMF models, in which the mass of the fictitious a meson is adjusted to nu- 
clear matter and ground-state properties of finite nuclei, whereas free values 
are used for the masses of the oo and p mesons. 

The point-coupling Lagrangian Eq. ([8]) does not include isovector-scalar 
terms. In the meson-exchange picture this channel is represented by the 
exchange of an effective 5 meson, and its inclusion introduces a proton- 
neutron effective mass splitting and enhances the isovector spin-orbit po- 
tential. Although the spin-orbit strength has a relatively well-defined value, 
the distribution between the scalar and vector channels is not determined 
by ground-state data. To reduce the number of adjustable parameters, the 
isovector-scalar channel may be omitted from an energy density functional 
that will primarily be used for the description of low-energy nuclear struc- 
ture. 

In general the strength parameters of the interaction terms in Eq. ([8|) 
are functions of the nucleon 4-current: 

f = ^^/^^ = pu^" , (9) 

where is the 4-velocity defined as (1 — v'^)~'^^'^{l,v). In the rest-frame 
of homogeneous nuclear matter: d = 0. The single-nucleon Dirac equation, 
the relativistic analogue of the Kohn-Sham equation, is obtained from the 
variation of the Lagrangian with respect to ^: 

[7^(ia^-S'^-S^)-(m + S5)]^ = 0, (10) 
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with the nucleon self-energies defined by the foUowing relations: 



T.'^ = av{p,)f + e^-^^A^^ (11) 
= as{Pv)Ps - Ss^Ps (13) 

TV 



= aTv{Pv)f ■ (14) 



In addition to the contributions of the isoscalar-vector four-fermion inter- 
action and the electromagnetic interaction, the isoscalar-vector self-energy 
includes the "rearrangement" terms S^, arising from the variation of 
the vertex functionals as, ay, and axv with respect to the nucleon fields 
in the density operator p. The inclusion of the rearrangement self-energy is 
essential for energy-momentum conservation and the thermodynamical con- 



sistency of the model [23|, l26|, [23]. E5 and Sj,^ denote the isoscalar-scalar 
and isovector-vector self-energies, respectively. 

In the relativistic density functional framework the nuclear ground state 
is represented by the self-consistent mean-field solution of the system of 
equations pOh - ()14p . with the isoscalar and isovector 4-currents and scalar 
density: 

N 



j/i = {4>o\i^iiii^\(l)o) = tpki/iipk , (15) 

k=i 

N 

J/i = (</'o|^7m^V'|0o> = tpklfiTtpk , (16) 

k=l 

N 

PS = (0o|^^|0o) = "^vltpkipk, (17) 

k=l 

where ipk are Dirac spinors, and the sum runs over occupied positive-energy 
single-nucleon orbitals, including the corresponding occupation factors v^. 
The single-nucleon Dirac equations are solved self-consistently in the "no- 
sea" approximation that omits the explicit contribution of negative-energy 
solutions of the relativistic equations to the densities and currents. Vac- 
uum polarization effects are implicitly included in the adjustable density- 
dependent parameters of the theory. 

To determine the density dependence of the coupling functionals as, 
ay, and arv one could start from a microscopic (relativistic) equation of 
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state (EoS) of symmetric and asymmetric nuclear matter, and map the cor- 
responding nucleon self-energies on the mean- field self-energies Eqs. (jlip - 
(jl4p that determine the single- nucleon Dirac equation (jlOp . This approach 
has been adopted, for instance, in RMF models based on Dirac-Brueckner- 



Hartree-Fock self-energies in nuclear matter [23|, |2J, |25[ , or on in- medium 



chir al p erturbation theory (ChPT) calculations of the nuclear matter EoS 



37], l38[. In general, however, energy density functionals determined directly 
from a microscopic EoS do not provide a very accurate description of data in 
finite nuclei. The reason, of course, is that a calculation of the nuclear mat- 
ter EoS involves approximation schemes and includes adjustable parameters 
that are not really constrained by nuclear structure data. The resulting 
bulk properties of infinite nuclear matter (saturation density, binding en- 
ergy, compression modulus, asymmetry energy) do not determine uniquely 
the parameters of nuclear energy density functionals, which usually must 
be further fine-tuned to ground-state data (masses and/or charge radii) of 
spherical nuclei. 

In a phenomenological construction of a relativistic energy density func- 
tional one starts from an assumed ansatz for the medium dependence of the 
mean- field nucleon self-energies, and adjusts the free parameters directly to 
ground-states data of finite nuclei. This procedure was used, for instance, 
in the construction of the relativistic density-dependent interactions TW- 



99 [2a], DD-MEl [23], DD-ME2 [28(, PKDD [40!], PKOl [4l|], DD-PCl [SJ]. 



3. Adjusting parameters to masses: the empirical functional DD- 
PCl 

In an analysis of relativistic nuclear dynamics [i^] , modern high-precision 
nucleon- nucleon (NN) potentials (Argonne Vis, Bonn A, CD-Bonn, Idaho, 
Nijmegen, Viow-fc) were mapped on a relativistic operator basis, and the 
corresponding relativistic nucleon self-energies in nuclear matter were cal- 
culated in Hartree-Fock approximation at tree level 0|. A very interesting 
result is that, at moderate nucleon densities relevant for nuclear structure 
calculations, all potentials yield very similar scalar and vector mean fields 
of several hundred MeV magnitude, in remarkable agreement with standard 
RMF phenomenology: at saturation density a large and attractive scalar 
field T,s ~ —400 MeV, and a repulsive vector field T,y « 350 MeV. The 
different treatment of short-distance dynamics in the various NN potentials 
leads to slightly more pronounced differences between the corresponding self- 
energies at higher nucleon densities. Generally, however, all potentials pre- 
dict a very similar density dependence of the scalar and vector self-energies. 
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In the chiral effective field theory framework, in particular, these self-energies 
are predominantly generated by contact terms that occur at next-to-leading 
order in the chiral expansion. 

Of course at the Hartree-Fock level these NN potentials do not yield 
saturation of nuclear matter. Nevertheless, the corresponding self-energies 
can be used as the starting point in the modeling of medium dependence of 
a relativistic nuclear energy density functional. The density functional DD- 
PCl which is representative of the class of semi-empirical REDFs and 
will be used in illustrative calculations throughout this work, was adjusted 
starting from the Hartree-Fock isoscalar scalar and vector self-energies of 
the Idaho N^LO potential [iS]. The strength and density dependence of 
the interaction terms of the Lagrangian Eq. ([8]) were parameterized by the 
following ansatz: 

asip) = as + {bs + csx)e-^s'' , (18) 
av{p) = av + ftye-"^^^, (19) 

where x = p/psat, and psat denotes the nucleon density at saturation in sym- 
metric nuclear matter. In the isovector channel the corresponding Hartree- 
Fock nucleon self-energies, obtained by directly mapping microscopic NN 
potentials on a relativistic operator basis, are presently not available. There- 
fore, as it was done in the case of functionals based on finite-range meson- 



exchange interactions TW-99 [26], DD-MEl [27|], DD-ME2 [28|, and PKOl 



41[, the density dependence of the isovector- vector coupling function was 
modeled based on the results of Dirac-Brueckner calculations of asymmetric 
nuclear matter 24i |: 

arvip) = brve-'^^^'' ■ (20) 

The parameters of a nuclear EDF can be constrained by the choice of 
the nuclear matter (symmetric and asymmetric) equation of state. The 
functional DD-PCl was eventually fine-tuned to experimental binding en- 
ergies of finite nuclei and, since the calculated nuclear masses are not very 
sensitive to the nuclear matter saturation density, this quantity, together 
with the compression modulus and the Dirac mass, were kept fixed. The 
saturation density psat = 0.152 fm^'^ is in accordance with values predicted 



by most modern relativistic functionals, such as DD-MEl 23], and DD- 
ME2 [28]. From these functionals also the Dirac effective nucleon mass was 
taken: mj, = m + S5 = 0.58m. The Dirac mass is closely related to the 
effective spin-orbit single-nucleon potential, and empirical energy spacings 
between spin-orbit partner states in finite nuclei determine a relatively nar- 
row interval of allowed values: 0.57 < m*j^/m < 0.61. In Ref. [sH] it was 
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shown that, to reproduce experimental excitation energies of isoscalar giant 
monopole resonances, point-couphng interactions require a nuclear matter 
compression modulus close to 230 MeV, and thus K^o = 230 MeV was used 
for DD-PCl. Nuclear structure data do not constrain the nuclear matter 
EoS at high nucleon densities. Therefore, in addition to psat-, i^*d^ ^^"^ ^oo, 
two additional points on the E{p) curve in symmetric matter were fixed to 
the microscopic EoS of Akmal, Pandharipande and Ravenhall fi^j, based on 
the Argonne Vis NN potential and the UIX three-nucleon interaction. This 
EoS has extensively been used in studies of high-density nucleon matter and 
neutron stars. 

The isovector channel of the energy density functional determines the 
density dependence of the nuclear matter symmetry energy 

S2{p) = 04 + -X"(/0 - Psat) + -1^^° (P - Psatf H • (21) 



sat 



The parameter po characterizes the linear density dependence of the symme- 
try energy, and AKq is the isovector correction to the compression modulus. 
Experimental masses, unfortunately, do not place very strict constraints on 



the parameters of the expansion of 52 (p) 46| , but self-consistent mean- field 
calculations show that binding energies can restrict the values of 52 at nu- 
cleon densities somewhat below saturation density, i.e. at p w 0.1 fm~^. 
Additional information on the symmetry energy can be obtained from data 
on neutron skin thickness and excitation energies of giant dipole resonances. 
Although values of neutron radii are available only for a small number of 
nuclei and the corresponding uncertainties are large, recent studies have 
shown that relativistic effective interactions with volume asymmetry 04 in 
the range 31 MeV < 04 < 35 MeV predict values for neutron skin thickness 
that are consistent with data, and reproduce experimental excitation ener- 



gies of isovector giant dipole resonances (cf . Ref. [47| and references therein) . 
For DD-PCl, therefore, the volume asymmetry was fixed at 04 = 33 MeV, 
and the symmetry energy was varied at a density that corresponds to an 
average nucleon density in finite nuclei: (p) = 0.12 fm~^. 

Until recently the standard procedure of fine-tuning non-relativistic or 
relativistic nuclear density functionals, was to perform a least-squares adjust- 
ment of a given set of free parameters simultaneously to a favorite nuclear 
matter EoS and to ground-state properties of about ten to twelve spheri- 
cal closed-shell nuclei. Deformed systems have generally not been included 
in the adjustment of parameters because calculation of deformed nuclei is 
computationally much more intensive, especially in a multi-parameter fit. 
Ground-state data of closed-shell nuclei, however, include long-range cor- 
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relations that cannot be absorbed into global functionals, e.g. correlations 
determined by the coupling to low-energy collective vibrations. The exci- 
tation energies and the structure of low-lying collective vibrational modes 
crucially depend on the details of single-nucleon levels in the vicinity of 
the Fermi surface and, therefore, the corresponding ground-state correla- 
tions obviously cannot display a smooth dependence on nucleon number. 
On the other hand, it is well known that energy density functionals or, at 
the level of practical application, self-consistent mean-field models provide 
a much better description of deformed, open-shell nuclei. The reason is that 
the mean-field approach includes the mechanism of spontaneous symme- 
try breaking that generates the most important ground-state correlations in 
deformed nuclei: quadrupole correlations in the p/i-channel and monopole- 
correlations in the pp-channel [3, S] ■ Additional rotational corrections that 
arise from the restoration of rotational symmetry vary rather smoothly over 
large mass intervals, and therefore can be included implicitly in the den- 
sity functional. The same reasoning applies to corrections originating from 
the restoration of particle number in the regime of strong pairing realized 
in well-deformed nuclei. It also appears that that those features of an ef- 
fective interaction that determine surface properties are better constrained 
by binding energies of deformed systems, as compared to spherical nuclei 
[49[. The functional DD-PCl was directly adjusted to binding energies of 
axially symmetric deformed nuclei in the mass regions A ~ 150 — 180 and 
A ~ 230 — 250. Similar approaches, employing data on both spherical and 
deformed nuclei, have recently been adopted in th e op timization of nuclear 
energy density functionals of Skyrme type 0, 0, m- 

Calculated masses of finite nuclei are primarily sensitive to the three 
leading terms in the empirical mass formula: volume, surface and symmetry 
energy 

B.E. = a^A + a,^2/3 ^ q / > + . . . ^ (22) 

where a^, and 04 correspond to the volume binding energy, surface en- 
ergy, and symmetry energy, respectively, at saturation density in nuclear 
matter. One can, therefore, generate families of effective interactions that 
are characterized by different values of a„ , and 04 (or symmetry energy at 
a lower density, as explained above), and determine which parametrization 
minimizes the deviation from the empirical binding energies. Of course, 
if a functional is adjusted by varying the volume, symmetry, and surface 
energies, the parameters that determine these quantities will generally be 
correlated because of Eq. (j22p . When only a small number of nuclei is con- 
sidered, satisfactory results can be obtained with various, in general linearly 
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dependent combinations of parameters. The parameters of the functional 
DD-PCl were thus determined in a careful comparison of the predicted 
binding energies with data, for a set of 64 nuclei with A ~ 150 — 180 and 
A w 230 - 250. 

The ground states of 64 axially deformed nuclei were calculated in the 
self-consistent mean field approximation. Pairing correlations were treated 
in the BCS constant-gap approximation with empirical pairing gaps (5-point 
formula), and the pairing model space included two major oscillator shells 
{2hLOo) above the Fermi surface. This approximation is justified because 
pairing correlations contribute only a very small portion to the total binding 
energy. In nuclei there is a clear separation of scales between the bulk 
contributions to the binding energies of the order of hundreds to more than 
thousand MeV, and the pairing energy of the order of ten MeV. To take into 
account pairing correlations in a calculation of the binding of nuclei close to 
/3-stability, therefore, it is sufficient to consider only the monopole part of the 
effective pairing interaction adjusted to experimental paring gaps. This is, 
of course, no longer true in studies of phenomena determined by structure 
effects in the vicinity of the Fermi surface, such as nuclear excitations or 
fission barriers, or in nuclei far from stability, where detailed properties of 
the effective interaction in the pairing channel become important. In this 
review we will also present examples with more realistic pairing interactions, 
e.g. a zero-range force or a separable version of the Gogny force. 

A careful analysis of deviations between calculated and experimental 
masses (mass residuals), showed a pronounced isospin and mass dependence 
of the residuals on the nuclear matter volume energy at saturation. To re- 
duce the absolute mass residuals to less than 1 MeV, and to contain their 
mass and isotopic dependence, strict constraints on the value of must be 
met. The narrow window of allowed values of the volume energy cannot be 
determined microscopically already at the nuclear matter level, but rather 
results from a fine-tuning of the parameters of the energy density functional 
to experimental masses. Calculated binding energies and charge radii are 
also sensitive to the choice of the surface coefficient that determines the 
surface energy and surface thickness of semi-infinite nuclear matter. For the 
optimal density functional DD-PCl, characterized by the following prop- 
erties at the saturation point: nucleon density psat = 0.152 fm~^, volume 
energy ay = —16.06 MeV, surface energy = 17.498 MeV, symmetry en- 
ergy 04 = 33 MeV, and the nuclear matter compression modulus Knm = 230 
MeV, in Fig. [1] we display the absolute deviations of the calculated binding 
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Figure 1: Absolute deviations of tlie calculated binding energies from the experimental 
values of the 64 axially deformed nuclei, as functions of the asymmetry coefficient (upper 
panel), and mass number (lower panel). Lines connect nuclei that belong to the isotopic 
chains shown in the legend. The theoretical binding energies are calculated using the 
relativistic density functional DD-PCl. 
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energies from the experimental values as functions of the isospin asymmetry 

2 {N - Zf 

and nucleon number. Positive deviations correspond to under-bound nuclei. 
The functional DD-PCl corresponds to the lowest value in the multi- 
parameter fit, and does not display any visible isotopic or mass dependence 
of the deviations of calculated masses. The absolute errors for all 64 axially 
deformed nuclei in the mass regions yl « 150 — 180 and A ~ 230 — 250 are 
smaller than 1 MeV. With stronger binding in symmetric nuclear matter 
(i.e. by increasing the absolute value of a^), the corresponding deviations of 
calculated binding energies become larger, and they also acquire a definite 
isotopic dependence. Reducing the absolute value of reverses the isotopic 
trend of the errors. 

The density functional DD-PCl has been further tested in calculations 
of properties of spherical and deformed medium-heavy and heavy nuclei, in- 
cluding binding energies, charge radii, deformation parameters, neutron skin 
thickness, and excitation energies of giant multipole resonances. Results 
have been compared with available data, and with predictions of the most 
successful finite-range meson-exchange relativistic effective interactions. In 
general, a very good agreement with data has been obtained except, per- 
haps, for the effect of overbinding of spherical closed-shell nuclei. DD-PCl, 
like virtually all relativistic mean-field models, is characterized by a rela- 
tively low effective nucleon mass and, when adjusted to masses of deformed 
nuclei, it overbinds spherical closed-shell systems. The well known problem 
of "arches" of mass residuals between shell closures could, in principle, be 
addressed by a functional that goes beyond the static mean-field approxima- 
tion and includes an explicit energy dependence of the nucleon self-energies. 
Very good results have been obtained for the excitation energies of giant 
monopole and dipole resonances in spherical nuclei, calculated using the rel- 
ativistic quasiparticle random-phase approximation based on the DD-PCl 
functional. The agreement with data validates the choice of the nuclear 
matter compressibility and symmetry energy for DD-PCl. The total num- 
ber of parameters is 10, similar to most non-relativistic Skyrme-type density 
functionals. The effective Lagrangian of DD-PCl contains only four inter- 
action terms except, of course, the Coulomb term (cf. Eq. (I8|)), and the 
10 parameters determine the density dependence of the strength functionals 
and reflect the complex nuclear many-body dynamics. 

In the following sections we will also employ DD-PCl in a series of spec- 
troscopic calculations utilizing models that go beyond the "mean-field" level. 
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and explicitly include collective correlations related to restoration of broken 
symmetries and fluctuations in quadrupole coordinates. Strictly speaking, 
those correlations that are treated explicitly should not be included in the 
density functional in an implicit way, that is, in the parameters adjusted 
to data that already include correlations. In future studies the solution 
could be to adjust the functional to pseudodata, obtained by subtracting 
correlation effects from experimental masses and, eventually, radii. In de- 
formed nuclei the dominant contribution to ground-state correlations is the 
rotational energy correction 52|, which is relatively simple to calculate. Ap- 
proximate methods have been developed that enable a systematic evaluation 
of correlation energies for the nuclear mass table [sS]. Starting from a set 
of pseudodata, one expects that the corresponding modifications of the pa- 
rameters of the energy density functional will be relatively small. However, 
even a small change in the relative contribution of various interaction terms 
could be the decisive factor in specific cases of soft potential energy surfaces, 
coexistence of prolate and oblate shapes, level ordering, etc. 



4. 3D Relativistic Hartree-Bogoliubov model with a separable pair- 
ing interaction 

Relativistic energy density functionals have been employed in analyses 
of properties of ground and excited states in spherical and deformed nuclei. 
For a quantitative analysis of open-shell nuclei it is necessary to consider 
also pairing correlations. Pairing has often been taken into account in a 
very phenomenological way in the BCS model with the monopole pairing 
force, adjusted to the experimental odd-even mass differences. In many 
cases, however, this approach presents only a poor approximation. The 
physics of weakly-bound nuclei, in particular, necessitates a unified and self- 
consistent treatment of mean-field and pairing correlations. This has led 
to the formulation and development of the relativistic Hartree-Bogoliubov 
(RHB) model [s^], which represents a relativistic extension of the conven- 
tional Hartree-Fock-Bogoliubov framework. The RHB model provides a uni- 
fied description of particle-hole (ph) and particle-particle {pp) correlations 
on a mean-field level by using two average potentials: the self-consistent 
mean field that encloses all the long range ph correlations, and a pairing 
field A which sums up the pp-correlations. The ground state of a nucleus is 
described by a generalized Slater determinant that represents the vac- 
uum with respect to independent quasiparticles. The quasiparticle operators 
are defined by the unitary Bogoliubov transformation of the single-nucleon 
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creation and annihilation operators: 



»t = T.^ikct + VikCi, (23) 
I 

where U and V are the Hartree-Bogohubov wave functions determined by 
the solution of the RHB equation. In coordinate representation: 

ho-m-X A \ ( Uk{r) \ _ ( Ukir) \ 

-A* -h}, + m + X )\ Vk{r) J ' \ Vk{r) J ' 

In the relativistic case the self-consistent mean-field corresponds to the 
single- nucleon Dirac Hamiltonian hu oi Eq. (|10p . m is the nucleon mass, 
and the chemical potential A is determined by the particle number subsidiary 
condition such that the expectation value of the particle number operator in 
the ground state equals the number of nucleons. The pairing field A reads 

^ab{r, r') = ^Yl ^-bcdir, r')AJcd(r-, r'). (25) 

c,d 

where Kibcd(^;'"') are the matrix elements of the two-body pairing interac- 
tion, and the indices a, b, c and d denote the quantum numbers that specify 
the Dirac indices of the spinor. The column vectors denote the quasiparticle 
wave functions, and E/^ are the quasiparticle energies. The dimension of the 
RHB matrix equation is two times the dimension of the corresponding Dirac 
equation. For each eigenvector {Uk,Vk) with positive quasiparticle energy 
Ek > 0, there exists an eigenvector {V^, U^) with quasiparticle energy — -E^. 
Since the baryon quasiparticle operators satisfy fermion commutation rela- 
tions, the levels Ej. and —Ej. cannot be occupied simultaneously. For the 
solution that corresponds to a ground state of a nucleus with even particle 
number, one usually chooses the eigenvectors with positive eigenvalues E^. 

The single-particle density and the pairing tensor, constructed from the 
quasi-particle wave functions 



p,d{ry) = Y,y:k{r)Vdk{r'), (26) 

fc>0 

ncd{ry) = Y.U*^{r)Vdk{r'), (27) 



A:>0 



are calculated in the no-sea approximation (denoted by A; > 0): the summa- 
tion runs over all quasiparticle states k with positive quasiparticle energies 
Ek > 0, but omits states that originate from the Dirac sea. The latter are 
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characterized by quasiparticle energies larger than the Dirac gap (~ 1200 
MeV). 

Pairing correlations in nuclei are restricted to an energy window of a few 
MeV around the Fermi level, and their scale is well separated from the scale 
of binding energies, which are in the range of several hundred to thousand 
MeV. There is no empirical evidence for any relativistic effect in the nuclear 
pairing field A and, therefore, a hybrid RHB model with a non-relativistic 
pairing interaction can be employed. For a general two-body interaction, 
the matrix elements of the relativistic pairing field read 

Aaipi,a2P2 = 2 ^ {aiPl,a2P2\V^^\a3P3,a4PA)a K.a:iP3,aiP4 , (28) 

where the indices {pi,P2,P3,P4 = f,g) refer to the large and small compo- 
nents of the quasiparticle Dirac spinors: 

U{r,s,t) = ( ) Vir, s,t) = ( ^^^f^'^'l ) . (29) 

In practical applications of the RHB model to finite open-shell nuclei, only 
the large components of the spinors Uk{r) and Vfc(r) are used to build the 
non-relativistic pairing tensor k in Eq. (I26p . The resulting pairing field reads 

(30) 



The other components: Ajg, A^j, and Aqq can be safely omitted [551 ]. 

In most applications of the RHB model [l^] the pairing part of the Gogny 
force [5fi] has been employed in the particle-particle {pp) channel: 

VPP(1,2) = ^ e-«''i-'"2)/^»)' (W^i + BiP'' - HiP^ - MiP^'P^) , (31) 

1=1,2 



with the set DIS 57| for the parameters Hi, Wi, Bi, Hi, and Mj (i = 1, 2). A 
basic advantage of the Gogny force is the finite range, which automatically 
guarantees a proper cut-off in momentum space. However, the resulting 
pairing field is non-local and the solution of the corresponding Dirac-Hartree- 
Bogoliubov integro-differential equations can be time-consuming, especially 
in the case 3D calculations for nuclei with triaxial shapes. 

In a series of recent articles 58, 5^, 60, 61 1 a separable form of the pairing 



force has been introduced for RHB calculations in spherical and deformed 
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nuclei. The force is separable in momentum space, and is completely de- 
termined by two parameters that are adjusted to reproduce in symmetric 
nuclear matter the bell-shape curve of the pairing gap of the Gogny force. 
The gap equation in the ^Sq channel reads 



and the pairing force is separable in momentum space: 

(A;| \k') = -Gp{k)p{k') . (33) 

By assuming a simple Gaussian ansatz p{k) = , the two parameters 

G and a have been adjusted to reproduce the density dependence of the 
gap at the Fermi surface, calculated with a Gogny force. For the DIS pa- 
rameterization 



57l | of the Gogny force the following values were determined: 
G = -728 MeVfm^ and a = 0.644 fm. When the pairing force Eq. (I33p is 
transformed from momentum to coordinate space, it takes the form: 

V{ru r2,r[,r'^) = G5{R- R') P(r)P(r')^ (1 - , (34) 

where R = ^ (?"i + T'2) and r = ri — r2 denote the center-of-mass and the 
relative coordinates, and P{r) is the Fourier transform of p{k): 

P{r) = ^e-^V4a^ . (35) 

(47ra2)3/2 

The pairing force has finite range and, because of the presence of the factor 
S {R — R'), it preserves translational invariance. Even though 6 {R — R') 
implies that this force is not completely separable in coordinate space, the 
corresponding anti-symmetrized pp matrix elements 

{aP\ V \-/d)^ = {aP\ V \-fd) - {a^\ V \h) , (36) 

can be represented as a sum of a finite number of separable terms in the 
basis of a 3D harmonic oscillator. 

The Dirac-Hartree-Bogoliubov equations ()24p are solved by expanding 
the nucleon spinors in the basis of a 3D harmonic oscillator in Cartesian 



coordinates |6l[ |. In this way both axial and triaxial nuclear shapes can be 
described. To obtain complete convergence, for medium heavy nuclei the 
basis must include at least NJ^""^ = 14 major oscillator shells, and for very 
heavy nuclei NJ^""^ = 16 major oscillator shells are necessary. The map of 



22 



the energy surface as a function of the quadrupole deformation is obtained 
by imposing constraints on the axial and triaxial quadrupole moments. The 
method of quadratic constraints uses an unrestricted variation of the func- 
tion 

{H)+ C2;.((Q2m>-<Z2m)', (37) 

M=0,2 

where (H) is the total energy, and {Q2fi) denotes the expectation value of 
the mass quadrupole operators: 

Q20 = - - and Q22 = - ■ (38) 

q2^ is the constrained value of the multipole moment, and 6*2^ the corre- 
sponding stiffness constant [sl]. For a self-consistent solution the quadratic 
constraint adds an extra force ^1^=0 2 ^tJ-Q2fi to the system, where = 
2C2^(((52^) — (l2^l)■ Such a force is necessary to keep the system at a point 
different from a stationary point. In general, the values of the quadrupole 
moments {Q2^) of the self-consistent solution coincide with the constrained 
values q2^ only at the stationary points. Moreover, the difference between 
the quadrupole moment (02^t) and the constrained value ^2/^ depends on 
the value of the stiffness constant, that is, smaller values of C2fi lead to 
larger deviations of the quadrupole moment from the corresponding con- 
strained value. Increasing the value of the stiffness constant, however, often 
destroys the convergence of the self-consistent procedure. This defici ency 
can be resolved by implementing the augmented Lagrangian method [63t], 
and this approach has been used in all constrained calculations presented in 
this work. 

Figs. [2] and [3] display the results of self-consistent 3D RHB calculations 
for the isotopic chains of Th, U, Pu, Cm, Cf, Fm, and No. The absolute 
deviations of the calculated binding energies from data [gS] (Fig. [2|) show 
an excellent agreement between theory and experiment. Except for the 
isotopes of No, all the calculated nuclei are slightly underbound, but the 
absolute differences between calculated and experimental binding energies 
is less than 1 MeV in all cases. The calculated ground-state quadrupole 
Q20J and hexadecapole Q^q, moments are compared with available data 
in Fig. [3l We notice that the values predicted by the DD-PCl functional 
reproduce in detail the isotopic trend of the empirical moments in the Th, U, 
Pu and Cm sequences, and are in very good agreement with the quadrupole 
moments of the Cf isotopes. No data on ground-state moments are available 
for the isotopes of Fm and No. The calculated 3D RHB quadrupole moments 



23 




are compared with values estimated in axially-symmetric calculations with 
the Lublin-Strasburg drop (LSD) model [g^ . 

The structure of the nucleus ^^'^Pu and its double-humped fission barrier 
has become a standard benchmark for models based on the self-consistent 
mean-field approach and the corresponding effective interactions or density 
functionals. In Fig.[l]we display the RHB triaxial quadrupole binding energy 
map of ^^'^Pu in the /3 — 7 plane (0 < 7 < 60*^), calculated with the DD- 



PCl energy density functional plus the pairing interaction Eq. (I34p [67]. 
The calculation has been carried out on a mesh of quadrupole deformation 
parameters with A/3 = 0.05 and A7 = 6^. All energies are normalized with 
respect to the binding energy of the absolute minimum, and the color code 
refers to the energy of each point on the surface relative to the minimum. 
Since the present implementation of the model does not include reflection 
asymmetric shapes, the potential energy surface (PES) is calculated only 
up to /3 < 1.3. For larger deformations, i.e. in the region of the second 
barrier, octupole deformations should be taken into account. The absolute 
minimum is calculated at /? = 0.28, 7 = O", and a second (super-deformed) 
valley is predicted around /? ~ 0.9. The axially symmetric barrier at /3 ~ 0.5 
is bypassed through the triaxial region, bringing the height of the barrier 
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Figure 3: Self-consistent RHB ground-state axial quadrupole and hexadecapole moments 
in comparison with data (open symbols), for the isotopic chains of Th, U, Pu, Cm, 
Cf, Fm, and No. For Fm and No the calculated quadrupole moments are compared with 
values predicted by the LSD model [g^ . 



much closer to the empirical value. This is shown more clearly in Fig. [5l 
where we plot the deformation energy curves and the inner barrier of ^^"^Pu 
as functions of the axial deformation (3. The two curves correspond to the 
axially-symmetric RHB calculation (solid), and to the projection on the /3- 
axis of the triaxial PES (dashed). The experimental values for the ground- 
state deformation, the barrier height, and the energy of the second minimum 
are taken from Refs. 6^, 6^, 7^, 7l| • One might notice a very good agreement 
between theory and available data. In particular, the inclusion of triaxial 
shapes lowers the inner barrier by ~ 2 MeV. Similar results have also been 
obtained in constrained self-consistent mean-field calculations using Skyrme 
functionals [2], and in the HFB-I-Gogny analysis of the actinide region |72j] it 
was shown that the inner barriers of the actinides are systematically lowered 
by up to 4 MeV when calculations included triaxial shapes. 

The 3D relativistic Hartree-Bogoliubov model, with the functional DD- 
PCl in the particle- hole channel and a separable pairing force in the particle- 
particle channel, enables very efficient constrained self-consistent triaxial 
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Figure 4: Self-consistent RHB triaxial quadrupole binding energy maps of Pu in the 
^-7 plane (0 < 7 < 60"). All energies are normalized with respect to the binding energy 
of the absolute minimum. The color code refers to the energy of each point on the surface 
relative to the minimum. 

calculations of binding energy maps as functions of quadrupole deformation 
in the /3 — 7 plane. The resulting single-quasiparticle energies and wave 
functions can be employed as microscopic input for the generator coordi- 
nate method configuration mixing of angular- momentum projected triaxial 
wave functions, or can be used to determine the parameters of the collec- 
tive Hamiltonian for vibrations and rotations: the mass parameters, the 
moments of inertia, and the collective potential. The solution of the cor- 
responding eigenvalue problem yields the excitation spectra and collective 
wave functions that are used in the calculation of electromagnetic transition 
probabilities. This approach will be illustrated in the next two sections. 
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Figure 5: (Color online) Deformation energy curves and the inner barrier of ■^*''Pu as 
functions of the axial deformation /?. The two curves correspond to the axially-symmetric 
RHB calculation (solid), and to the projection on the /3-axis of the triaxial PES (dashed), 
calculated with the functional DD-PCl. The experimental values for the ground-state 
deformation, the barrier height, and the energy of the second minimum are indicated, 
respectively, with an arrow, a symbol with error bars and three lines indicating the value 
and its errors. The data are taken from Refs. [6^. |6^. [toI. ItH . 

5. Beyond the mean-field approximation: restoring broken sym- 
metries and configuration mixing calculations 

Nuclear structure far from stability has become a subject of extensive 
experimental and theoretical studies. The variation of ground-state shapes 
in an isotopic chain, for instance, is governed by the evolution of shell struc- 
ture. Far from the /3-stability line, in particular, the energy spacings be- 
tween single-particle levels change considerably with the number of neutrons 
and/or protons. This can result in reduced spherical shell gaps, modifica- 
tions of shell structure, and in some cases spherical magic numbers may 
disappear. The reduction of a spherical shell closure is associated with the 
occurrence of deformed ground states and, in a number of cases, with the 
phenomenon of shape coexistence. 

A quantitative description of structure phenomena related to shell evolu- 
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tion necessitates the inclusion of many-body correlations beyond the mean- 
field approximation. The starting point is usually a constrained Hartree- 
Fock plus BCS (HFBCS), or Hartree-Fock-Bogoliubov (HFB) calculation 
of the potential energy surface with the mass quadrupole components as 
constrained quantities. When based on microscopic EDFs or effective inter- 
actions, such calculations comprise short- and long-range many-body corre- 
lations, and result in static symmetry-breaking product many-body states. 
Static mean-field models provide a description of bulk properties, such as 
masses and radii. To calculate energy spectra and transition probabilities, 
however, it is essential to include collective correlations that arise from sym- 
metry restoration and fluctuations around the mean-field minimum. Because 
they are susceptible to shell effects and vary with nucleon number, collective 
correlations cannot be be incorporated in a universal EDF. In deformed nu- 
clei, for instance, the rotational energy correction, that is, the energy gained 
by the restoration of rotational symmetry, is proportional to the quadrupole 
deformation of the symmetry-breaking state and can reach several MeV 
for a well deformed configuration. Fluctuations of quadrupole deformation 
also contribute to the correlation energy. Both types of correlations can be 
included simultaneously by mixing angular-momentum projected states cor- 
responding to different quadrupole moments. The most effective approach 
for configuration mixing calculations is the generator coordinate method 
(GCM), with multipole moments used as coordinates that generate the in- 
trinsic wave functions. 

In recent years several accurate and efficient models and algorithms, 
based on microscopic density functionals or effective interactions, have been 
developed that perform the restoration of symmetries broken by the static 
nuclear mean field, and take into account quadrupole fluctuations. Many in- 
teresting phenomena related to shell evolution have been investigated by em- 
ploying the angular-momentum projected GCM with the axial quadrupole 
moment as the generating coordinate, and with intrinsic configurations cal- 
culated in the HFB model with the finite-range Gogny interaction 



75l . l76l . l77l |. Very recently this approach has been extended to include full 



triaxial angular-momentum and particle-number projection [78|]. Another 
sophisticated structure model that takes into account collective correlations 
is based on axially constrained HF-I-BCS calculations with Skyrme effective 
interactions in the particle-hole channel and a density-dependent contact 
force in the pairing channel 7^, M, M, [si, l83l |. Particle numbers and rota- 
tional symmetry are restored by projecting self-consistent mean- field wave 
functions on the correct numbers of neutrons and protons, and on angular 
momentum. Finally, a mixing of the projected wave functions corresponding 
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to different quadrupole moments is performed with a discretized version of 
the generator coordinate method. The latest extension of this model that in- 



corporates triaxial angular-momentum projection was reported in Ref. [84]. 

In a series of recent articles we have expanded the framework of relativis- 
tic energy density functional to include correlations related to the restora- 
tion of broken symmetries and to fluctuations of collective variables. A 
model has been developed that uses the GCM to perform configuration mix- 



ing of angular-momentum [15|, and also particle-number projected [16] rel- 
ativistic wave functions. The geometry was restricted to axially symmetric 
shapes, and the intrinsic wave functions were generated from solutions of the 
relativistic mean-field + Lipkin-Nogami BCS equations, with a constraint 
on the mass quadrupole moment. In the first application |[85j], the GCM 
based on relativistic EDFs was employed in a study of shape transitions in 
Nd isotopes. This approach has been further developed in Refs. 3, 19] by 



implementing a model that includes triaxial angular-momentum projection. 

Here we will illustrate the "beyond mean-field" extension of the rela- 
tivistic EDF approach by considering a rather simple example of a GCM 
configuration mixing of angular-momentum projected wave functions, gen- 
erated from the self-consistent solutions of axially symmetric constrained 
RMF+BCS equations for ^^^^Sm. 

5.1. GCM mixing of angular-momentum projected states 

The generator coordinate method (GCM) is based on the assumption 
that, starting from a set of intrinsic symmetry-breaking states \(p{q)) which 
depend on a collective coordinate q, one can build approximate eigenstates 



of the nuclear Hamiltonian 64] 



\^a) = ^faiq,)\HQj)) ■ (39) 
j 

Here the basis states \4>{q)) are Slater determinants of single- nucleon states 
generated by solving the constrained relativistic mean-field + BCS equations 
with the mass quadrupole moment as the generating coordinate q. The 
axially deformed mean-field breaks rotational symmetry, so that the basis 
states \(p{q)) are not eigenstates of the total angular momentum. To be able 
to compare theoretical results with data, it is necessary to construct states 
with good angular momentum 

l^f = E fi''i9j)PMK m)) , (40) 
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where P^jj^ denotes the angular momentum projection operator 

Pmk = / dnDi*^in)R{n) . (4i) 

The integration is performed over the three Euler angles a, /3, and 7. 
Dljj^{n) = e-*^°d5^^(/3)e-^^T is the Wigner function, and R{n) = e"*"-^- 
e~^^-^y e~'''^'^'' is the rotation operator. The weight functions fa^{qj) are 
determined from the variation: 

that is, by requiring that the expectation value of the energy is stationary 
with respect to an arbitrary variation Sf^^. This leads to the Hill-Wheeler 
equation 

E /"""(^j) (^-^(ft)! ^^MK I'PiQj)) Ei P-l,K mi))) = . (43) 

The restriction to axially symmetric configurations (J^ |^ (<?)) = 0) simplifies 
the problem considerably, because in this case the integrals over the Euler 
angles a and 7 can be performed analytically, and for an arbitrary multipole 
operator Qx^ one thus finds 

(<^fe)| Qxt^P^MK \4>{<li)) = (2J + l)'-^^SM-^.SKo Jo^' dPsinP d%^{l3) 

mi)\Qx^.e-'^'^\^{q,)) . (44) 

This expression vanishes for odd values of angular momentum J and, there- 
fore, all projected quantities are defined only for even values of J. The norm 
overlap kernel 

J^'{m,qj) = {HQ^)\PMK\<PiQJ)) = 

(2J + l)^d±^SMoSKo r dpsinp d4*(/3) 
^ Jo 

mi)\e-'^'ymj)) , (45) 
can be evaluated by employing the generalized Wick theorem: 

nim, qj; p) = {<j){qi)\ e''^'^ |<^(g,)) = ±^ det Mab{qi,qf, P) ■ (46) 
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The overlap matrix is defined: 



■hfabiQi, Qj;P) = Ua{qi)Rab{qi, Qj', P)ub{qj) + Va{qi)Rab{qi, qj', P)vb{qj) , (47) 
where u and v denote the BCS occupation probabihties, the matrix R reads 

Rabiqi, qj;f3) = J i^iir; qi)e-'f^^yMr; qj)dr , (48) 

and il>a{r; qi) denotes the self-consistent intrinsic single-nucleon spinor at the 
generating coordinate g^. The Hamiltonian kernel 

n-'iqi^qj) = {ct>{qi)\HPUm3)) = 

(2 J + 1) ' SmoSko dl3 sin /3 ^5 iP) 

{ct>{qi)\He-'^'y\cl>{qj)) , (49) 

can be calculated from the energy functional provided the intrinsic densities 
are replaced by transition one-body matrices when evaluating the expression 

h{qi, qj;P) = (</>(gi)| ^e"'^""^ \<t>{qj)) = J dr Stot{r; m, qj,P) . (50) 

The basis states |^(gj)) are not eigenstates of the proton and neutron 
number operators Z and N. The adjustment of the Fermi energies in a 

BCS calculation ensures only that the average value of the nucleon number 
operators corresponds to the actual number of nucleons. Consequently, the 
wave functions are generally not eigenstates of the nucleon number 

operators and, moreover, the average values of the nucleon number operators 

are not necessarily equal to the number of nucleons in a given nucleus. To 
restore the correct mean values of the nucleon numbers, one can modify the 
Hill- Wheeler equation by replacing h{qi,qj;f3) with 

h'{qi, qj\P) = h{qi, qj; (3) - Xp [z{qi, qj; /3) - zq] - A„ [n{qi, qj; (3) - no] , (51) 

where 

z{qi,qj;P) = {HQi)\ Ze-'^'" \<t>{<lj)) , (52) 

and 

n{qi, qr, /3) = {4>{qi)\ Ne-'^'y |0(g,)) . (53) 

\p(^n) is the proton (neutron) Fermi energy, while zq and uq denote the 
desired number of protons and neutrons, respectively. 
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The Hill- Wheeler equation 



= EiY,^\q,,qj)fi{q,) , (54) 

3 j 

presents a generalized eigenvalue problem, and thus the weight functions 
faiQi) are not orthogonal and cannot be interpreted as collective wave func- 
tions for the variable q. It is useful to re-express Eq. ()54p in terms of another 
set of functions, g^iqi), defined by 

3 

With this transformation the Hill- Wheeler equation defines an ordinary 
eigenvalue problem 

= Eo.gi{qi) , (56) 

3 

with 

'H\q:,qj) = Y.^M'y^'\q,,qk)'H\qk.qi){M')-^/\quq^) . (57) 
k,l 

The functions gaiQi) ^ire orthonormal and play the role of collective wave 
functions. 

Once the weight functions fail) are known, it is straightforward to cal- 
culate all physical observables, such as transition probabilities and spec- 
troscopic quadrupole moments. The reduced transition probability for a 
transition between an initial state {Ji,ai), and a final state {Jf,af), reads 



B{E2;Jiai Jjaj) 



e2 



2Ji + l 



2 

Jf* 



qf,qi 



(58) 



5.2. Angular momentum projection and configuration mixing: ^^^5m 

As an illustration of GCM configuration mixing of angular-momentum 
projected states, we consider the case of the prolate deformed, axially- 
symmetric ground state of ^^^Sm. The intrinsic symmetry-breaking wave 
functions that are used in the configuration mixing calculation are obtained 
as solutions of the self-consistent relativistic mean- field equations, subject to 
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Figure 6: The RMF+BCS binding energy curve of ^^*Sm (dashed), and the corresponding 
angular-momentum projected {J'" = 0^,2+, 4+ , and 6^) energy curves, as functions of 
the axial deformation /3. 

constraint on the mass quadrupole moment. The interaction in the particle- 
hole channel is determined by the density functional DD-PCl, and a density- 
independent (5-force is used as the effective interaction in the particle-particle 
channel, supplemented with a smooth cut-off determined by a Fermi function 
in the single-particle energies [33]. Pairing correlations are treated within 
the BCS framework, and the pairing contribution to the total energy is given 
by 

= / ir)dr = {r)dr , (59) 

for protons and neutrons, respectively. Kp(n)(^) denotes the local part of the 
pairing tensor, and Vp = —321 MeV fm'^ and Vn = —308 MeV fm'^ are the 
pairing strength parameters adjusted to empirical pairing gaps. We note 
that while the BCS constant-gap approximation was employed in the fit of 
the parameters of the functional DD-PCl to experimental binding energies 
(cf. Sec.[3|), obviously constant gaps cannot be used in a constrained calcula- 
tion. Constant gaps that are adjusted to empirical values correspond only to 
the ground state and, therefore, cannot be used to calculate configurations 
that generally represent excited states. This is the reason for employing a 
zero-range pairing force in this example of angular momentum projection 
and configuration mixing calculation. 

The GCM basis is constructed from the self-consistent solution of con- 



33 



strained single-nucleon Dirac equations on a regular mesh in the generating 
coordinate - the mass quadrupole moment. The corresponding axial de- 
formation parameter spans the interval from (3 = —0.68 to /3 = 1, with a 
spacing A/3 = 0.04. The GCM basis thus consists of 43 intrinsic states. The 
large and small components of the Dirac spinors are expanded in a basis of 
axially symmetric oscillator eigenfunctions for = 14 major shells. 

The resulting binding energy curve of ^^^Sm is shown in Fig. [6l as func- 
tion of the axial deformation /?. The RMF+BCS mean-field calculation 
yields a prolate ground-state minimum at /3 = 0.32, more than 4.5 MeV 
deeper than the wide local minimum on the oblate side. One might notice 
that for prolate shapes the potential is rather stiff with respect to /3 around 
the minimum, and this means that the structures built on the ground state 
will be localized in this minimum, that is, there will be very little mixing with 
oblate structures for states with relatively low angular momenta. Fig. [6] also 
displays the corresponding angular-momentum projected {J'^ = 0"'", 2"*", 4"*", 
and 6"*") energy curves. At this stage we do not consider configuration mixing 
yet, and the projected energy of the \4>{q)) state reads 



are plotted in Fig. [71 The spherical configuration is a pure 0+ state (7^'^=°(0, 
0) = 1), and the maxima of the projected norm overlap kernels for higher 
angular momenta are correspondingly shifted to larger deformations. 

Coming back to the projected energy curves in Fig. O we notice that, 
since the spherical configuration is already a pure 0"*" state, there is no energy 
gain for J'^ = O"*" at /? = 0. The rotational energy correction at the prolate 
minimum /3 = 0.32 is more than 2.5 MeV, and even larger in the region 
of the oblate minimum. The spherical point /3 = is not included in the 
plots of E-^ (q) for J > 2. Namely, for J / the quantities 'H"^(0, 0) and 
AA'^(0, 0) are so small, that their ratio in Eq. (I60p cannot be determined 
accurately. For higher values of the angular momentum several additional 
configurations close to the spherical point are also characterized by very 
small values of the projected norm overlap kernel. These configurations can 
safely be omitted from the projected energy curves, because on the one hand 
the angular momentum projection becomes inaccurate at these points, and 
on the other hand the corresponding angular momentum projected states 
would not play any role in configuration mixing calculations. 




(60) 



where the projected norm overlap kernels 




(61) 
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Figure 7: Projected norm overlap kernels Af''{q,q) as a functions of the axial quadrupole 
deformation /3. 



The next step in the solution of the Hill- Wheeler equation ()54p is the 
diagonalization of the norm overlap kernel Af'^{qi,qj) 

^^M-\qi,qj)uk{qj) = rikUkiqi) ■ (62) 



Since the basis functions \(i){qi)) are not linearly independent, many eigen- 
values Uk will be very close to zero. The corresponding eigenfunctions Uk{qi) 
are rapidly oscillating and carry very little physical information. However, 
due to numerical uncertainties, their contribution to 'H'^{qi,qj) can be large, 
and these states should be removed from the basis. From the remaining 
states one builds the collective Hamiltonian 

nif = -^^Y.uMn'{q^,q,)ui{q,) , (63) 

which is subsequently diagonalized 

Y,ni[gf^ = Eigi'^ . (64) 
k,l 
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Figure 8: Angular-momentum projected GCM results for the excitation energies and 
B(E2) values (in Weisskopf units) of the lowest two bands in ^^''Sm, in comparison with 
available data. 



The solution determines both the ground state energy, and the energies of 
excited states, for each value of the angular momentum J. The collective 
wave functions ga{q), and the weight functions fi{q), are calculated from 
the norm overlap eigenfunctions 

/ 

and 

fi{q^) = Y,^Mq^)■ (66) 

The GCM excitation energies and the corresponding B(E2) values for 
the two lowest bands in ^^^Sm: the ground-state band and the /3-band, are 
displayed in Fig. [HI The results of the angular-momentum projected (AMP) 
configuration mixing calculation are compared with available data. Since 
the present calculation is restricted to axial symmetry, it cannot describe 
structures based on the 7 degree of freedom. Considering that, for a given 
EDF and the effective pairing interaction, the AMP+GCM calculation is 
parameter- free, the agreement with data is remarkable. One might notice 
the excellent agreement of the transition probabilities, calculated with the 
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Figure 9: Amplitudes of the angular-momentum projected GCM collective wave functions 
\gi{P)?' {J"' = 0+,2+,4+, and 6+) for the ground-state band (left) and the /3-band (right) 
in ^s^Sm. 



bare proton charge, with the experimental values for transitions within the 
ground-state band, and the transition from the band-head of the /3-band to 
the ground-state band. The calculated excitation energy of the /3 band-head 
is above the experimental 0^ , indicating that the potential is probably too 
stiff with respect to /3, and the corresponding moment of inertia appears 
to be considerably lower than the effective empirical value. The difference 
between the calculated and empirical moments of inertia is much less pro- 
nounced in the ground-state band. 

In Fig. [9]we plot the amplitudes of the collective wave functions 
Eq. (|65p for the two lowest GCM states of each angular momentum J'^ = 
0^,2"^,4~^, and 6^. The amplitudes in the left panel correspond to the 
states of the ground-state band, whereas the collective functions of the /3- 
band are plotted in the right panel. For both bands the collective functions 
are localized in the prolate well, and only for the state 0^ the wave function 
displays a small amount of oblate admixture. 

The treatment of collective correlations in REDF-based structure models 
has been extended in Ref. [16] to axially symmetric GCM configuration mix- 
ing of angular-momentum and particle-number projected states, and very re- 
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cently a model has been developed that includes triaxial angular-momentum 



projection 18|, ll9|. The latter does not involve particle number projection. 
It has to be emphasized that, while GCM configuration mixing of axially 
symmetric states has been implemented by several groups and routinely used 
in nuclear structure studies, the application of this method to triaxial shapes 
presents a much more involved and technically difficult problem. Only the 
most recent advances in parallel computing and modeling have enabled the 



implementation of models [8J, |78| , based on triaxial symmetry-breaking in- 
trinsic states, that are projected on particle number and angular momentum, 
and finally mixed by the generator coordinate method. This implementa- 
tion is equivalent to a seven-dimensional GCM calculation, mixing all five 
degrees of freedom of the quadrupole operator and the gauge angles for pro- 
tons and neutrons. The numerical realization, however, is very complex, 
and applications to medium-heavy and heavy nuclei are still computation- 
ally too demanding and time-consuming. This is true even for a model that 



does not include particle number projection [19|]. In addition, the use of 
general EDFs with an arbitrary dependence on nucleon densities in GCM- 
type calculations, often leads to discontinuities or even divergences of the 



energy kernels as a function of deformation 86|, |87|. Only for a specific type 
of density dependence a regularization method can be implemented 8^, 89], 
that corrects the energy kernels and removes the discontinuities and diver- 
gences. In the next section we therefore review an approximation to the full 
GCM approach, based on the microscopic REDF framework, that includes 
rotational symmetry restoration and takes into account triaxial quadrupole 
fiuctuations in arbitrary heavy nuclei. 

6. Collective Hamiltonian in five dimensions based on relativistic 
EDFs 

In an alternative approach to five-dimensional quadrupole dynamics that 
restores rotational symmetry and allows for fiuctuations around the triax- 
ial mean-field minima, a collective Bohr Hamiltonian can be formulated, 
with deformation-dependent parameters determined by microscopic self- 
consistent mean-field calculations. There are two principal approaches to 
derive the collective Hamiltonian starting from a microscopic framework 
based on an effective inter-nucleon interaction or energy density functional: 
(i) the adiabatic approximation to the time-dependent HFB theory (AT- 
DHFB) 90|, and (ii) the generator coordinate method (GCM) with the 



Gaussian overlap approximation (GOA) |9ll. l92l. l93l|. With the assumption 



that the GCM overlap kernels can be approximated by Gaussian functions 
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[63], the local expansion of the kernels up to second order in the non-locality 
transforms the GCM Hill- Wheeler equation into a second-order differential 
equation - the Schrodinger equation for the collective Hamiltonian. The 
kinetic part of this Hamiltonian contains an inertia tensor [9^], and the po- 
tential energy is determined by the diagonal elements of the Hamiltonian 
kernel, and also includes zero-point energy (ZPE) corrections {s^. 

The dynamics of the collective Bohr Hamiltonian is determined by the 
vibrational inertial functions and the moments of inertia 96|]. For these 
quantities either the GCM-GOA (Yoccoz masses [o^]) or the ATDHFB ex- 
pressions (Thouless-Valatin masses [o^) can be used. The Thouless-Valatin 
masses have the advantage that they also include the time-odd components 
of the mean-field potential and, in this sense, the full dynamics of a nuclear 
system. In the GCM approach these components can only be included if, 
in addition to the coordinates gj, the corresponding canonically conjugate 
momenta pi are also taken into account, but this is obviously a very compli- 
cated task. In many applications a further simplification is thus introduced 
in terms of cranking formulas |99|, |95|], that represent the perturbative limit 
for the Thouless-Valatin masses, and the corresponding expressions for ZPE 
corrections. This approximation was applied in recent studies using models 
based both on t he G ogny interaction jinnl . and Skyrme energy den- 



sity functionals |102l . l93l |. The approximate inclusion of Thouless-Valatin 



corrections to the mass parameters and mome nts o f inertia of the Skyrme- 
based Boh r Hamiltonian was discussed in Ref. 102l |. In a recent systematic 
study [103l | of low-energy nuclear structure at normal deformation, based on 
the Hartree-Fock-Bogoliubov theory extended by the generator coordinate 
method and mapped onto a five-dimensional collective quadrupole Hamilto- 
nian, the Thouless-Valatin moments of inertia were used, whereas the crank- 
ing approximation was used for the quadrupole mass parameters. Using the 
Gogny DIS interaction, even-even nuclei with proton numbers Z = 10 to 
Z = 110 and neutron numbers N < 200 were calculated. 

Here we review a recent implementation for the solution of the eigenvalue 
problem of a five-dimensional collective Hamiltonian for quadrupole vibra- 
tional and rotational degrees of freedom, with parameters determined by 
constrained sel f-con sist ent relativistic Hartree-Bogoliubov calculations for 
triaxial shapes [loi, [13, [loi, [lO^, [lO^, Ifi? | . 



6.1. Collective Hamiltonian 

The general Bohr collective model for the description of quadrupole 
collective states, including a detailed discussion of the model's kinemat- 
ics, has recently been reviewed in Ref. [93|]. Nuclear excitations deter- 
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mined by quadrupole vibrational and rotational degrees of freedom can be 
treated simultaneously by considering five quadrupole collective coordinates 
a^, fj, = —2,-1,..., 2 that describe the surface of a deformed nucleus: 
R = i?o[l + '^f^o^^J.Y2f2■ To separate rotational and vibrational motion, 
these coordinates are usually parameterized in terms of two deformation 
parameters /3 and 7, and three Euler angles {(j), 9, ip) = Q that define the 
orientation of the intrinsic principal axes in the laboratory frame 

= Dl,{n)(3cos^ + i= [D%in) + Dl_,in)] /3sin7 , (67) 

where -D^^ is the Wigner function. The three terms of the classical collective 
Hamiltonian, expressed in terms of the intrinsic variables /?, 7 and Euler 
angles 

^cou = rvib(/3,7) +7;ot(/3,7,f^) + Vcoii(/3,7) , (68) 
denote the contributions from the vibrational kinetic energy: 

Tvib = \BfSf3$^ + pB^^P^ + ^f^^B^^^^ , (69) 



the rotational kinetic energy: 

3 

2 



1 ^ 

%ot = ::Y.^kUjl (70) 



k=l 



and the collective potential energy Vcoii(/3,7)- The mass parameters 

-B-y^y, and the moments of inertia Ik depend on the quadrupole defor- 
mation variables /3 and 7. 

After quantization the classical kinetic energy of the Hamiltonian Eq. ()68p : 

reads: 

The kinetic energy tensor takes the block diagonal form: 
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with the vibrational part of the tensor 



vib 



B 



(74) 



In general, the rotational part is a complicated function of the Eulcr angles 
but, using the quasi-coordinates related to the components of the angular 
momentum in the body-fixed frame, it takes a simple diagonal form 



(-Brot)ifc — ^iA, k — 1,2,3 , 

with the moments of inertia expressed as 

Xfc = ABkP^ sin2(7 - 2A;7r/3) . 



(75) 



(76) 



This particular functional form is motivated by the fact that all three mo- 
ments of inertia vanish for the spherical configuration {/3 = 0) and, addi- 
tionally, Iz and ly vanish for axially symmetric prolate (7 = 0*^) and oblate 
(7 = 60°) configurations, respectively. The resulting determinant reads 



detB = detSvib ■ detSrot = 4wr/3^ sin^ 87 , 



(77) 



where w = B^j^B^^ — B^^ and r = BiB2B^. The quantized collective 
Hamiltonian can hence be written in the form: 



H = Tvib + Trot + KoU 



(78) 



with 

Tvib = - 

+ 

and 



2^ 



d 



d d 



dp dpy w 



d 



1 



P sin 87 



d 



d 



1 ^ P 



k=i 



(79) 



(80) 



where Jk denotes the components of the angular momentum in the body- 
fixed frame of a nucleus. V^oU is the collective potential. The Hamiltonian 
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describes quadrupole vibrations, rotations, and the coupling of these collec- 
tive modes. The determinant Eq. (j77p determines the volume element in 
the collective space: 



/f POO pZtt n 

dTcoii= J dndToVwr = J J ^7] sin 87! J dn^/IiPF , (81) 

and the quantized Hamiltonian Eq. (|78p is hermitian with respect to the 
collective measure Eq. ()8ip . 

The eigenvalue problem of the general collective Hamiltonian Eq. ()78p 
can be solved by a direct numerical solution of a system of partial differential 
equations using finite-difference methods, or by employing an expansion of 
eigenfunctions in terms of a truncated basis in the collective Hilbert space. 
The basis functions depend on the deformation variables /3 and 7, and the 
Euler angles (j), 6 and ^|J. In the latter case the eigenvalue problem reduces 
to a simple matrix diagonalization, and the main task is the construction 
of an appropriate basis for each value of the angular momentum quantum 
number. In the implementation developed in Ref. 171] we employed the basis 
expansion approach. 

The diagonalization of the Hamiltonian yields the excitation energies 
and collective wave functions: 

^I'f^(/3,7,fi)= i^iMlXKi^)- (82) 

K&AJ 

The angular part corresponds to a linear combination of Wigner functions 



and the summation in Eq. (|82p is over the allowed set of the K values: 

^j^(02 J for ^mod2 = 

\ 2, 4, . . . , J - 1 for J mod 2 = 1. ^ ^ 

Using the collective wave functions Eq. ()82p . various observables can be cal- 
culated and compared with experimental results. For instance, the quadrupole 
E2 reduced transition probability: 

B{E2; aJ ^ a'J') = ^^\{a'J'\\Mim\»J)W (85) 

where Ai{E2) is the electric quadrupole operator. For the 7W(ii^2) matrix 
elements the current implementation of the model uses a local expression 
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in the collective deformation variables 108|. This approximation is justified 



in the case of large overlaps between different vibrational amplitudes jlOOl ]. 
but may be less suited for transitions between states with a rather small 
overlap, e.g. for transitions between super-deformed bands and bands at 
normal deformation. 

The shape of a nucleus can be characterized in a qualitative way by the 
average values of the invariants /3^, cos 87, as well as their combinations. 
For example, the average value of the invariant in the state |aJ): 

{I3^)j^ = {K\li''\0 = E / /?'Kx(/3,7)l'rfro , 

and the average values of the deformation parameters /3 and 7 in the state 
I a J) are calculated from: 

(/3)ja = VW)7a, (87) 

1 (/?3cos37)ja 
(7)ja = o arccos ■ 



3 ^Wh^mTe. 

The mixing of different intrinsic configurations in the state \al) can be deter- 
mined from the distribution of the projection K of the angular momentum 
/ on the z axis in the body-fixed frame: 

Nk = & / |V;^,i^(/3,7)|'/3^|sin37|6?/3d7, (89) 



where the components V'o_ft-(/5,7) are defined in Eq. ([82]) . For large de- 
formations the K quantum number is to a good approximation conserved. 
Consequently, only one of the integrals Eq. (|89p will give a value close to 1. 
A broader distribution of Nk values in the state \aJ) provides a measure of 
mixing of intrinsic configurations. 

6.2. Microscopic parameters of the collective Hamiltonian 

The entire dynamics of the collective Hamiltonian is governed by the 
seven functions of the intrinsic deformations /3 and 7: the collective poten- 
tial, the three mass parameters: Bpi^, Bp^, B^^, and the three moments 
of inertia X^. These functions are determined by the choice of a partic- 
ular microscopic nuclear energy density functional or effective interaction. 
The entire map of the energy surface as function of the quadrupole defor- 
mation is obtained by imposing constraints on the axial and triaxial mass 
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quadrupole moments, as described in Sec [H The quasiparticle wave func- 
tions and energies, generated from constrained self-consistent solutions of 
the RHB model, provide the microscopic input for the parameters of the 
collective Hamiltonian. 

In the simplest approximation the moments of inertia are calculated from 
the Inglis-Belyaev formula: 



E 



(90) 



where k denotes the axis of rotation, the summation runs over proton and 
neutron quasiparticle states = /3|/3j|<I>), and \^) represents the quasi- 
particle vacuum. The mass parameters associated with the two quadrupole 
collective coordinates qo = {Q20) and (72 = {Q22) are calculated in the crank- 
ing approximation: 



where 



M(n),fiA<l0,q2) = J2 



-^ri)-^(3)-^(i) 



mQ2^.\ij){^j\Q2u\^) 



{Ei + Ej 



(91) 



(92) 



The collective energy surface includes the energy of the zero-point mo- 
tion, and this quantity has to be subtracted. The collective zero-point en- 
ergy (ZPE) corresponds to a superposition of zero-point motion of individual 
nucleons in the single-nucleon potential. In the general case, the ZPE cor- 
rections on the potential energy surfaces depend on the deformation. The 
ZPE includes terms originating from the vibrational and rotational kinetic 
energy, and a contribution of potential energy 



^^(^0,92) = AKib(9o,92) + AV;ot(go,92) + AVpot(go,92) 



(93) 



The latter is much s mall er than the contribution of the kinetic energy, and 
is usually neglected [100]. Simple prescriptions for the calculation of vibra- 
tional and rotational ZPE have been derived in Ref. 95|]. Both corrections 
are calculated in the cranking approximation, i.e. on the same level of 
approximation as the mass parameters and the moments of inertia. The 
vibrational ZPE is given by the expression: 



AKib(9o,92) = jTr 



^(3)^(2) 



(94) 
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The rotational ZPE is a sum of three terms: 



AVrot{qo,q2) = AV-2-2{qo,q2) + AV-i-i{qo,q2) + AVu{qo,q2), (95) 
with 

AV(go,g2) = -r-rz 7 r • (96 

4 A^(3),Mi.(90,g2) 

The individual terms are calculated from Eqs. (|96p and (|92p . with the 
intrinsic components of the quadrupole operator defined by: 

Q21 = -2iyz , Q2-1 = -2xz , Q2-2 = 2ixy . (97) 



The potential V^oU in the collective Hamiltonian Eq. (178p is obtained by 
subtracting the ZPE corrections from the total mean-field energy: 

Koii(go, 92) = -E^tot(go, 92) - AKib(9o, 92) - AV;ot(go, ^2) • (98) 



As an example in Fig. [TU] we display the spectrum of collective states 
of ^'''^Pu. Starting from constrained self-consistent solutions of the RHB 
equations, i.e. using single-quasiparticle energies and wave functions that 
correspond to each point on the energy surface shown in Fig. the pa- 
rameters of the collective Hamiltonian are calculated as functions of the 
deformations (3 and 7. The diagonalization of the Hamiltonian yields the 
excitation spectrum shown in Fig. [TOl shown in comparison to data for the 
three lowest positive-parity bands at normal deformation, and the lowest 
vr = + super-deformed band of ^^'^Pu. In addition to the yrast ground-state 
band, in deformed and transitional nuclei excited states are also assigned 
to (quasi) /3 and 7 bands. This is done according to the distribution of the 
projection K of the angular momentum / on the z axis of the body-fixed 
frame Eq. (j89p . Excited states with predominant K = 2 components in the 
wave function are assigned to the 7-band, whereas the /3-band comprises 
the states above the yrast characterized by dominant K = components. 
K = states are assigned to the super-deformed band based on the calcu- 
lated average value of the deformation parameter /3 Eq. ([87)) . 

The Inglis-Belyaev (IB) moments of inertia Eq. (I90p of the collective 
Hamiltonian have been multiplied by a common factor so that the calcu- 
lated energy of the 2^ state coincides with the experimental value. This 
scale parameter reflects the well known fact that the IB expression predicts 
effective moments of inertia that are smaller than empirical values. In the 
calculation of the spectrum of ^^'^Pu we have thus followed the prescription 



45 




Figure 10; The low-energy spectrum of Pu calculated with the DD-PCl relativistic 
density functional (left), compared with data (right) for the three lowest positive-parity 
bands at normal deformation, and the lowest vr = + superdeformed band. 



of Ref. lOOd where, by comparing the more reahstic, but also computa- 
tionally more involved, Thouless-Valatin (TV) moments of inertia with the 
IB values as functions of the axial deformation for superdeformed bands in 
the A ~ 190 mass region, it was shown that the TV correction to the per- 
turbative IB expression is almost independent of deformation, and does not 
include significant new structures in the moments of inertia. It was thus sug- 
gested that the moments of inertia to be used in the collective Hamiltonian 
can be simply related to the IB values through the minimal prescription: 
^k{Q) = ^k^{Q){^ + ct)) where q denotes the generic deformation parameter, 
and a is a constant that can be determined in comparison to data. In the 
present case a = 0.32 for ^'^'^Pu. 

When the IB effective moment of inertia is renormalized to the empirical 
value, the excitation spectrum of the collective Hamiltonian determined by 
the functional DD-PCl is in very good agreement with the available data for 
the ground-state band, (3 and 7 bands, and even the lowest super-deformed 
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band SD-1. Compared to the corresponding experimental sequence, the 
position of the 7 band is predicted at somewhat lower excitation energy, 
and this might indicate that the theoretical PES is probably too soft in 7. 
The /3-band is calculated at slightly higher energy compared to experiment, 
and the predicted position of SD-1 is within the experimental error bounds. 
Very few data are available on electromagnetic transition rates in ^^"Pu. In 
fact, except for the lifetime of the 2f state, only the lifetimes of K-isomers 
have been measured but these include configurations not contained in our 
collective model space. Therefore, in Fig. [10] we only display the calculated 
B(E2) values, in Weisskopf units (W.u.), for the transition 2^ — t- 0^ and 
from the band-heads of the /3 and 7 bands to the ground-state band. It 
should be emphasized that besides the renormalization of the moment of 
inertia, the calculation is completely parameter-free, i.e. by using structure 
models based on self-consistent mean-field single-particle solutions, physical 
observables, such as transition probabilities and spectroscopic quadrupole 
moments, are calculated in the full configuration space and there is no need 
for effective charges. Using the bare value of the proton charge in the elec- 
tric quadrupole operator Ai{E2), the transition probabilities between eigen- 
states of the collective Hamiltonian can be directly compared to data. 

6.3. Illustrative calculation: evolution of triaxial shapes in Pt isotopes 

Most deformed nuclei display axially-symmetric prolate ground-state 
shapes, but some regions of the nuclide chart are characterized by the oc- 
currence of oblate deformed and triaxial shapes. One of the examples is the 
A ~ 190 mass region, where both prolate to oblate shape transitions, as well 
as triaxial ground-state shapes have been predicted. An extensive analysis 
of this region has re cently bee n performed using non-relativistic Skyrme and 



Gogny interactions 109l . Ill0l |. The self-consistent Hartree-Fock-Bogoliubov 
model has been used to study the evolution of the ground-state shapes of 
Yb, Hf, W, Os and Pt isotopes. In particular, it has been shown that the 
isotopic chains with larger Z-numbers in this mass region display a tendency 
toward triaxial shapes. Here we present the 3D RHB binding energy maps 
for the sequence of even-A Pt isotopes with neutron numbers in the interval 
from = 108 to A^ = 126, calculated with the DD-PCl energy density 
functional plus the pairing interaction Eq. (134p . 

In Tab. [1] we list the calculated values of the (3 and 7 deformation pa- 
rameters for the absolute minima of the potential energy surfaces (PES). 
One can follow the transition from the prolate deformed ^^^Pt, through the 
region of triaxially deformed i88-i98pt isotopes, to the slightly oblate ^"^Pt, 
and finally the spherical 202-204p|. iggtopes. The ground-state /3-deformation 
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13 



Figure 11: Self-consistent RHB triaxial quadrupole binding-energy maps of the even-even 
isotopes i90-200pt in the /3 - 7 plane (0 < 7 < 60°). All energies are normalized with 
respect to the binding energy of the absolute minimum (red dot) . 
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A 


186 


188 


190 


192 


194 


196 


198 


200 


202 


204 




0.30 


0.28 


0.19 


0.18 


0.15 


0.13 


0.12 


0.08 








7 


0° 


9° 


34° 


34° 


34° 


31° 


33° 


60° 


0° 


0° 



Table 1: Calculated values of the /3 and 7 deformation parameters for the absolute minima 
of the potential-energy surfaces (PES) of even-A Pt isotopes with 186 < A < 204. 



steadily decreases as the number of neutrons increases and approaches the 
closed-shell at = 126. In order to analyze the nature of shape transi- 
tion in the Pt isotopic chain, in Fig. [TT] we display the self-consistent RHB 
quadrupole binding energy maps of the even-A i90-200p^ isotopes in the 
/3 - 7 plane (0° < 7 < 60°). Ah energies are normalized with respect to the 
binding energy of the absolute minimum, and the color code refers to the 
energy of each point on the surface relative to the minimum. The PES of 
190-198 are 7-soft, with shallow minima at 7 ~ 30*^. The nucleus ^^''Pt 
displays a slightly oblate minimum, signaling the shell-closure at A = 126. 

As an illustrative example for the microscopic origin of the triaxial 
ground-state deformations, we consider the nucleus ^^^Pt. The formation of 
deformed minima can be related to the occurrence of gaps or regions of low 
single-particle level density around the Fermi surface. In Figs. [12] and [13] 
we plot the proton and neutron single-particle energy levels in the canonical 
basis for ^^^Pt. Solid curves correspond to levels with positive parity, and 
short-dashed curves denote levels with negative parity. The long-dashed 
(yellow) curve corresponds to the Fermi level. The leftmost and the right- 
most panels display prolate and oblate axially-symmetric single-particle lev- 
els, respectively, whereas the middle panel shows the single-particle levels as 
functions of 7 for the fixed value of the axial deformation = 0.18. This 



type of plot has been introduced in Ref. |llll |. and it enables the identifica- 
tion of K quantum numbers of triaxial sing l e-particle levels in the limits of 
axial symmetry at 7 = O'' and 7 = 60*^ 1091 110 . II2I ]. 



In Fig. [12] we notice the occurrence of a gap between the proton single- 
particle levels in the vicinity of the Fermi surface around 7 = 30'^. The 
energy gap predominantly results from the down-sloping of one particular 
single-particle orbital, originating from the spherical shell, as the de- 
formation parameter 7 increases from 7 = 0° to 7 = 60'^. This result is 
in agreement with the findings of Ref. \lO^. The corresponding neutron 
single-particle levels, shown in Fig. [TS] also display a region of low level 
density around the Fermi surface at 7 ~ 30", although the gap is somewhat 
less pronounced in comparison to the proton gap. 
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protons 

prolate p=0.175 oblate 




Figure 12: Proton canonical single-particle energy levels of Pt. Solid curves denote 
levels denote levels with positive parity and short-dashed curves levels with negative par- 
ity. The long-dashed (yellow) curve corresponds to the Fermi level. The leftmost and 
the rightmost panels display prolate and oblate axially symmetric single-particle levels, 
respectively. The middle panel shows the single-particle levels as functions of 7 for the 
fixed value of the axial deformation — 0.18. 
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Figure 13: Same as described in the caption to Fig. [12] but for neutron single-particle 
levels. 



The low-energy spectra of Pt and Pt, obtained by diagonalization 
of the collective Hamiltonian based on the DD-PCl energy density func- 
tional plus the pairing interaction Eq. (|34p . are displayed in Figs. O and 
[TSl The calculated ground-state bands and (quasi) 7-bands are compared 
with the corresponding sequences of experimental states. Both the theoret- 
ical excitation energies and B(E2) values are in very good agreement with 
data. In the case of Pt isotopes it was not necessary to renormalize the 
effective moments of inertia, i.e. the spectra shown in Figs. [14] and [15] do 
not include any additional scaling parameter. As already emphasized, tran- 
sition probabilities are calculated in the full configuration space with bare 
proton charges. In particular, one might notice the excellent result for the 
predicted excitation energy of the bandhead of the 7-band in both nuclei, 
as well as the very good agreement with the experimental B(E2) values for 
transitions between the 7-band and the yrast band. This result indicates 
that the DD-PCl potential has the correct stiffness with respect to the 7 
degree of freedom. A similar agreement with data is also obtained in the 
calculation of low-energy spectra for the other Pt isotopes whose PES are 
displayed in Fig. [TTJ 
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Figure 14: The low energy spectrum of ^^'^Pt calculated with the DD-PCl relativistic 
density functional (left) compared with the data (right) for the excitation energies and 
intraband and interband B(E2) values (in Weisskopf units). 



4.0 
3.5 
3.0 
2.5 
2.0 



> 

0) 



S 1.5 



1.0 
0.5 
0.0 




Figure 15: Same as described in the caption to Fig. 1141 but for the isotope ^'^''Pt. 
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The five-dimensional collective Hamiltonian for quadrupole vibrational 
and rotational degrees of freedom, with parameters determined by con- 
strained self-consistent relativistic mean-field calculations for triaxial shapes, 
has also been employed in studies of microscopic signatu res of gro und-state 
shape phase transitions in Nd isotopes around = 90 105 . 10d |. and Ba 



and Xe nuclei in the mass A ~ 130 region jl07f |. A detailed c ompa rison with 
avail able data and with predictions of the analytical X(5) [ll3 | and E(5) 



114( 1 models, has shown that the microscopic theoretical framework based 



on relativistic EDFs describes not only general features of shape transitions, 
but also singular properties of excitation spectra and transition rates at the 
points of first- and second-order quantum shape phase transitions. 

There are possible improvements and extensions of the model that has 
been reviewed in this section. For instance, because in most cases the Inglis- 
Belyaev formula yields effective moments of inertia that are lower than em- 
pirical values, all the calculated relative excitation energies had to be scaled 
with respect to the experimental energy of the 2f states. The moments of 
inertia can be improved by including the Thouless-Valatin dynamical rear- 
rangement contributions. For the rotational degrees of freedom for which 
the collective momenta are known, the inertia parameters can be obtained 
from the solutions of cranked RHB equations. For the deformation coordi- 
nates qo and q2 the situation is more complicated, because the corresponding 
momentum operators Pq and P2 have to be calculated from the solution of 
Thouless-Valatin equations [os'] at each deformation point. Because cranking 
breaks time-reversal symmetry, in both cases the inclusion of pairing cor- 
relations necessitates calculations in the full relativistic Hartree-Bogoliubov 
framework, including time-odd components of the nucleon self-energies. 



7. Summary and outlook 

A wealth of new data from radioactive-beam facilities, the exciting phe- 
nomenology of nuclear astrophysics, and recent theoretical developments 
in related fields, have prompted important advances in theoretical nuclear 
structure physics during the last decade. The objective of this field is to 
build a consistent microscopic framework that will, on the one hand, bridge 
the gap between the underlying theory of strong interactions and the phe- 
nomenology of finite nuclei and, on the other, provide a unified description of 
bulk properties, excitations and reactions across the entire chart of nuclides. 

Even though ab initio approaches, starting from a microscopic nuclear 
Hamiltonian that accurately reproduces scattering and few-body data, have 
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been very successful in the description of relatively light nuclei up to oxy- 
gen isotopes, and large-scale semi-microscopic shell model calculations are 
performed for medium-heavy and even some heavy nuclei in the vicinity of 
closed shells, at present the only comprehensive approach to nuclear struc- 
ture is provided by the framework of energy density functionals. The advan- 
tages of using EDFs in the description of structure phenomena are evident 
already at the basic level of implementation - the self-consistent mean-field 
method: an intuitive interpretation of mean-field results in terms of intrinsic 
shapes and single-particle states, calculations are performed in the full model 
space of occupied states (no distinction between core and valence nucleons, 
no need for effective charges), and the universality of EDFs that enables 
their applications to all nuclei throughout the periodic chart. The latter 
feature is especially important for extrapolations to regions of exotic short- 
lived nuclei far from stability for which few, if any, data are available. For 
spectroscopic applications, however, the EDF-based approach must be ex- 
tended beyond the static mean-field level and models must be developed that 
include collective correlations related to the restoration of broken symme- 
tries in finite nuclei, and take into account fluctuations of collective variables 
around mean-field minima. 

Relativistic energy density functionals (REDFs) have their origin in the 
highly successful relativistic mean-field (RMF) phenomenological models in- 
troduced by Walecka and Serot 0, 0] , and later applied and further devel- 
oped by many groups. More recently, this framework has been reinter- 
preted by analogy to relativistic Kohn-Sham density functional theory. It 
has been realized that the original meson-exchange forces used in RMF 
models present only one of the possible representations of the effective in- 
medium inter- nucleon interactions and, moreover, one that does not present 
any particular advantage at low energies characteristic for nuclear binding 
and low-lying excitations. Functionals have thus been developed that are 
expressed in terms of ground-state nucleon four-currents and scalar densities 
only, with short-distance correlations and long-range dynamics represented 
either by higher order powers of the currents and densities, or encoded in the 
medium (nucleon density) dependence of the coupling functions of interac- 
tion Lagrangians. The corresponding structure models have been applied to 
studies of a variety of phenomena in spherical and deformed nuclei, extend- 
ing over the whole mass table and to systems with extreme isospin values. 
The illustrative calculations presented in this work have been performed us- 
ing the relativistic energy density functional DD-PCl [34] for which, start- 
ing from microscopic nucleon self-energies in nuclear matter, the parameters 
were fine-tuned in a careful fit to experimental binding energies of 64 axially 
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deformed nuclei in the mass regions A Ri 150 — 180 and A ^ 230 — 250. For 
quantitative calculations in open-shell nuclei it also necessary to consider 
pairing correlations and, when used in the relativistic Hartrec-Bogoliubov 
(RHB) framework together with a pairing force separable in momentum 
space, the functional DD-PCl provides an excellent description of ground- 
state properties. The corresponding self-consistent QRPA calculations re- 
produce the excitation energies of giant multipole resonances. 

In this work we have also reviewed the latest extensions of the REDF 
framework that include the treatment of collective correlations. By restoring 
symmetries broken by the static mean-field and considering fluctuations of 
collective deformation variables, REDF-based models have been developed 
that can be employed in detailed spectroscopic studies, including predic- 
tions for excitation spectra and electromagnetic transitions. For axially 
deformed nuclei this approach has been illustrated by GCM configuration 
mixing calculations of angular-momentum projected relativistic mean-field 
wave functions. The GCM excitation energies and the corresponding B(E2) 
values for the two lowest bands in ^^^Sm have been discussed in comparison 
with available data. In an approximation to the full GCM treatment of the 
five-dimensional quadrupole dynamics, a collective Bohr Hamiltonian has 
been formulated, with deformation-dependent parameters determined by 
constrained microscopic self-consistent RHB calculations of triaxial energy 
surfaces in the /3 — 7 plane. The entire map of the energy surface as function 
of the quadrupole deformation is obtained by imposing constraints on the 
mass quadrupole moments. The quasiparticle wave functions and energies, 
generated from constrained self-consistent solutions of the RHB model, pro- 
vide the microscopic input for the parameters of the collective Hamiltonian: 
the collective potential, the three mass parameters: Bjsis, B^^, -B-,-,, and the 
three moments of inertia I^- The implementation of this complex, REDF- 
based, microscopic collective model has been exemplified in a study of the 
evolution of triaxial shapes in Ft isotopes. The DD-PCl based 3D RHB 
calculation provides a simple microscopic interpretation of the occurrence of 
triaxial shapes in terms of proton and neutron single-particle levels, and the 
calculated ground-state bands and (quasi) 7-bands are in very good agree- 
ment with the corresponding sequences of experimental states, both for the 
excitation energies and B(E2) values. 

In the remainder of this section we outline some of the most important 
challenges for the framework of (relativistic) density functionals. 

An important issue is the development of a series (ladder) of accurate and 
controlled approximation for the exchange-correlation terms of the energy 
density functional, analogous to the "Jacob's Ladder" of Coulomb Density 
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Functional Theory |115l |. One possible approach is to develop the exchange- 
correlation functional from first principles by incorporating known exact 
constraints, another is empirical and optimizes a parametric ansatz by ad- 
justing it to a set of data. In the context of nuclear density functionals this 
topic has recently been discussed in the review of Ref. [6(], and a possible 
non-empirical approach to climbing the rungs of the ladder of approxima- 
tions toward the universal functional has been indicated. 

A related topic is the link between the framework of nuclear EDFs and 
the underlying theory of strong interaction - low-energy QCD, which will 
include both nuclear matter and finite nuclei. At low energies characteristic 
for n uclear binding, QCD is realized as a theory of pions coupled to nucleons 



1161]. The basic concept of a low-energy effective field theory (EFT) is the 



separation of scales: the long-range physics (pion exchange) is treated ex- 
plicitly, whereas short-distance interactions, that cannot be resolved at low 
energy, are replaced by contact terms. In a non-empirical approach to nu- 
clear NEDFs, the EFT-based derivation of exchange-correlation functionals 
in principle allows for error estimates, and provides a power counting scheme 
that separates long- and short-distance dynamics. 

Even a non-empirical universal energy density functional will have to be 
fine-tuned to data on medium-heavy and heavy nuclei. This is because such 
a functional contains a number of low-energy constants that determine the 
strength of the leading short-range interactions. These constants cannot be 
adjusted already from scattering and few-body data to an accuracy that en- 
ables a quantitative description of structure phenomena in complex nuclear 
systems. In principle any complete set of low-energy data, for instance ex- 
perimental masses, can be used to fine-tune the EDF. However, when only 
a small number of nuclei is considered, satisfactory least-squares fits can 
be obtained with different, in general linearly dependent combinations of 
parameters. Moreover, ground states of spherical nuclei that have mostly 
been used to adjust functionals or effective interactions, include collective 
correlations that cannot be absorbed in a universal functional. Another im- 
portant issue concerns functionals that are used in models that go beyond 
the mean-field level and include collective correlations. For instance, if rota- 
tional energy corrections and quadrupole fluctuations are treated explicitly 
by angular momentum projection and configuration mixing, they should 
not at the same time implicitly be included in the functional, i.e. through 
parameters adjusted to data that already include these correlations. There- 
fore, the parameters of such functionals must be adjusted to pseudodata, 
obtained by subtracting correlation effects from experimental masses and 
radii. 
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In the development of EDF-based structure models that include collec- 
tive correlations through symmetry restorations and configuration mixing, 
relativistic functionals face the same challenges as their non-relativistic coun- 
terpart, i.e. Skyrme-type functionals. For these models to be able to make 
spectroscopic predictions in medium-heavy and heavy nuclei, often charac- 
terized by soft potential energy surfaces, it is important to build accurate 
and efficient algorithms that perform a complete restoration of symmetries 
broken by the static mean field (translational, rotational, particle number), 
and take into account fluctuations around the mean-field minima for very 
general shapes. 

An interesting recent development, which goes beyond the scope of the 
present review, is the extension of the r elativistic (Q)RPA to the quasipar- 
ticle time-blocking approximation 117 . 118 . 1191 . Il2d . 121]. This approach 



takes into account effects of particle-vibrational coupling and, therefore, 
enables a quantitative analysis of single-particle excitations in odd-mass nu- 
clei and vibrational excitations. In particular, the dynamics of particle- 
vibrational coupling leads to an increase of the level density near the Fermi 
surface, i.e. to an enhancement of the nucleon effective mass. The RQRPA 
extended by the coupling to collective vibrations generates spectra with a 
multitude of two-quasiparticle-plus-phonon states, that are important in the 
description of damping phenomena characteristic for giant multipole excita- 
tions, as well as in studies of low-energy modes in neutron-rich nuclei. 

Finally, interesting results could also be obtained by expanding nuclear 
energy density functionals to include non-nucleonic degrees of freedom. For 
instance, it has been shown that relativistic density functionals provide a 
natural framework for the descri ption of hy pernuclear single-particle spectra 
based on chiral SU(3) dynamics [l2i.[l23j. 
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